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Abstract 

In the study of the heat transfer in the Bohzniann theory, the basic problem is to 
construct sohitions to the following steady problem: 

v\/,F = ^Q{F,F), (x,«)er!xR3, (0.1) 
F{x,v)\n(a:)-v<o = t^B F {x ,v'){n{x) ■ v')dv' , xedQ, (0.2) 

Jn{x)-v'>0 

where i7 is a bounded domain in R'', 1 < d < 3, Kn is the Knudsen number and fie = 
2it9^(x) 6xp[— 2e(x) \ ^ Maxwellian with non-constant (non-isothermal) wall temperature 
0{x). Based on new constructive coercivity estimates for both steady and dynamic cases, 
for |0 — ^ol ^ (5 ^ 1 and any fixed value of Kn, we construct a unique solution Fs to 
()0.ip and l|0.2p . continuous away from the grazing set and exponentially asymptotically 
stable. This solution is a genuine non equilibrium stationary solution differing from a 
local equilibrium Maxwellian. As an application of our results we establish the expansion 
Fs = + '^'^1 + 0{5^) and we prove that, if the Fourier law holds, the temperature 
contribution associated to Fi must be linear, in the slab geometry. This contradicts 
available numerical simulations, leading to the prediction of breakdown of the Fourier law 
in the kinetic regime. 
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1 Introduction and notation 

According to the Boltzmann equation (jl.2p . a rarefied gas confined in a bounded domain, 
in contact witli a tliermal reservoir modeled by ()0.2p at a given constant temperature 6 
(isothermal), has an equilibrium state described by the Maxwellian 

constant x e"'"'''/^^, 

and it is well known [35l[10l[18l[12l[39l[22l|38] that such an equilibrium is reached exponentially 
fast, at least if the initial state is close to the above Maxwellian in a suitable norm. 

If the temperature 9 at the boundary is not uniform in (j0.2p . such a statement is not true 
and the existence of stationary solutions and the rate of convergence require a much more 
delicate analysis, because rather complex phenomena are involved. For example, suppose that 
the domain is just a slab between two parallel plates at fixed temperatures 9- and 9+ with 
9^ < 9-^-. Then, one expects that a stationary solution is reached, where there is a steady flow 
of heat from the hotter plate to the colder one. The approach to the stationary solution may 
involve convective motions, oscillations and possibly more complicated phenomena. Even the 
description of the stationary solution is not obvious, and the relation between the heat flux 
and the temperature gradient, (e.g. the Fourier law (jl.ip ). is not a priori known. A first 
answer to such questions can be given confining the analysis to the small Knudsen number 
regime (i.e., Kn — )■ in (10. ip ). where the particles undergo a large number of collision per unit 
time and a hydrodynamic regime is established. In this case, it can be formally shown, using 
expansion techniques [U [H] , that the lowest order in Kn is a Maxwellian local equilibrium 
and the evolution is ruled by macroscopic equations, such as the Navier-Stokes equations. In 
particular, the heat flux vector q turns out to be proportional to the gradient of temperature, 
as predicted by the Fourier law 

q=-K{9)VJ (1.1) 

with the heat conductivity n{9) depending on the interaction potential. This was first obtained 
by Maxwell and Boltzmann \28\ [7] which relates the macroscopic heat flow to the microscopic 
potential of interaction between the molecules. The rigorous proof of such a statement was 
given in [16^ [T7] in the case of the slab geometry and provided that 9^ — 9- is sufficiently 
small (uniformly in the Knudsen number Kn). This is a special case of a problem which has 
received recently a large attention in the Statistical Mechanics community, the derivation 
of the Fourier law from the microscopic deterministic evolutions ruled by the Newton or 
Schrodinger equation or from stochastic models [8l |30l El H] . 

The aim of this paper is to analyze the thermal conduction phenomena in the kinetic 
regime. This problem was studied in the slab geometry, for small Knudsen numbers in \16\ ITT] 
and for large Knudsen number (Kn — )• cxd) in [30]. Here we are interested in a general domain 
and in a regime where the Knudsen number Kn is neither small nor large. In this regime, the 
only construction of solutions to (jO.ip and (j0.2p we are aware of, was achieved in [3] in a slab 
for large ^+ — 9-, with techniques closer in the spirit to the DiPerna-Lions renormalized 
solutions |13t I14j (see also [38] and references quoted therein). However, the uniqueness and 
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stability of such solutions are unknown. To develop the quantitative analysis we have in 
mind, the theory of the solutions close to the equilibrium ([351 IMl ttOl [22]) is better suited. 
For this reason we confine ourselves to temperature profiles at the boundary which do not 
oscillate too much. More precisely, we will assume that the temperature 6{x) on dil. is given 
by 9{x) = 6q + 5'd{x) with 6 a small parameter and 'd{x) a prescribed bounded function on 
(90 such that 

sup < 1. 

xedn 

This will allow us to use perturbation arguments in the neighborhood of the equilibrium at 
the uniform temperature ^o- 

We shall consider the Boltzmann equation 

dtF + vVF=^Q{F,F), (1.2) 

with F{t, X, v) the probability density that a particle of the gas at time t is in a small cell of 
the phase space x centered at {x,v). Here is a bounded domain in R*^, d = 1, 2, 3 with 
a smooth boundary dQ. The function F is required to be a positive function on $7 x R^ such 
that //f^xR3 Fdxdv is fixed for any t. The right hand side of (jl.2p . Q{F, G), is the Boltzmann 
collision operator (non-symmetric) 

Q{F,G) = [ dv^ [ dujB{v-v^,uj)F{v'^)Giv') 



du* / dujB{v — v*,u>)F{v*)G{v) 

= Qgain(i^,G)-gioss(i^,G), (1.3) 

where B{v,uj) = |upgo • with < 7 < 1 (hard potential), < qQ{uj ■ j^) < C\uj ■ |^| 

(angular cutoff) is the collision cross section and v', are the incoming velocities in a binary 
elastic collision with outgoing velocities v, and impact parameter w: 



V 



' = V — uj[{v — v^,) ■ u], vl = + u)[{v — v^,) ■ u]. (1.4) 



The contact of the gas with thermal reservoirs is described by suitable boundary con- 
ditions. We confine ourselves to the simplest interesting case of the diffuse refiection (10. 2p . 
although more general boundary data could be studied [TO]. On dQ (supposed to be a C^- 
smooth surface with external normal n{x) well defined in each point x E di}) we assume the 
condition: 

F{t,x,v) = fig{x,v) / F{t,x,v'){n{x) ■ v'}dv' , (1.5) 

Jn{x)-v'>0 

for X G and n{x) ■ v < 0, where fie is the Maxwellian at temperature 6, 



/^^(^>-) = ^;^exp 



29{x) 



(1.6) 



normalized so that 

fj,Q{x,v){n{x) ■ v}dv = 1. (1.7) 

n{x)-v>0 
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Throughout this paper, is a connected and bounded domain in R'^, for d = 1, 2, 3 and 
the velocity v G R*^ and -0 G R^~'^ such that 



V = {vi,--- ,Vd,Vd+i,--- ,V3) = {v,v). 



(1.8) 



We denote the phase boundary in the phase space i7 x R^ as 7 = dO, x R'^, and spht it into 
the outgoing boundary 7_|_, the incoming boundary 7_, and the grazing boundary 70: 



7+ = {{x,v) e dn xh^ 

7_ = {{x,v) £ dnxK^ 
70 = {{x,v) € dn xK^ 



n{x) ■ V > 0}, 
n(x) • V < 0}, 
n{x) ■ V = 0}. 



The backward exit time tY,{x,v) is defined for {x,v) G $7 x R'^ 

th{x, v) = inf{ t >0 : X — tv e 90}, 



(1.9) 



and X]3{x, v) = x — t\3{x, v)v G d^l. Furthermore, we define the singular grazing boundary 7q , 
a subset of 70, as: 



7^ = {(x,u)G7o : tb{x,-v)^0 and tbix,v)j^O}, 
and the discontinuity set in Q x R'^: 

S = 70 U G n X R^ : (xb(x, f ), f ) G 7^ }. 

We will use the short notation fig for the Maxwellian 
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27r[^o + 6^x)]' 



■ exp 



2[eo + 5^{x)] 



(1.10) 



(1.11) 



(1.12) 



Moreover, to denote the global Maxwellian at temperature 9o, HOo-, we will simply use the 
symbol ji: 

Since (|1.7|) is valid for all 5, we have 



/ 



li{v){n{x) ■ v}dv = 1. 



(1.13) 



i(a::)-D>0 

We denote by L the standard linearized Boltzmann operator 

Lf = -^[Q{j^,VJ^f) + QiVJ^f,i^)] = Hv)f - Kf 

= v{v)f - I li{v,v^)f{v^)dv^, 

yR3 



(1.14) 



with the collision frequency i'{v) = J /j^a^sa B{v — v^,uj)^{v^)dojdv^ ~ {1 + for < 7 < 1. 
Moreover, we set 



r(/l, /2) = -^^(^//^/l, VA^/2) = rgam(/l, /2) " ^IossUiJt)- 



(1.15) 



Finally, we define 

P^f{x,v) = sj li{v) I f{x,v')y^fi{v'){n{x) ■ v')dv'. (1-16) 

Jn(x)-v'>0 

Thanks to (jl.lSp . P-yf, viewed as function on {v G R'^ | v ■ n{x) > 0} for any fixed x E dO., 
is a L^-projection with respect to the measure |n(x) • -ul for any boundary function / defined 
on 7+. 

We denote || • ||oo either the L°°(ri x R'^)— norm or the L°°(r2)— norm in the bulk, while 
I • loo is either the L°°{dQ. x R^)— norm or the L°°(c?r2)— norm at the boundary. Also we 
adopt the Vinogradov notation: X <Y \s equivalent to \X\ < CY where C is a constant 
not depending on X and Y. We subscript this to denote dependence on parameters, thus 
X <aY means \X\ < CaY. Denote (v) = ^l + \v\^. Our main results are as follows. 

Theorem 1.1 There exists 6o > such that for < 6 < 6o in and for all M > 0, 

there exists a non-negative solution Fs = Mfi + y/JIfs > with //f^xR^ fsy/Jj'dxdv = to the 
steady problem jlO. 1\} and W.^) such that for all ^ ( < j:^, /3 > 4, 

\\{vfe<\^\'fs\\oo + \{vfe'^^-\'fs\o.<S. (1.17) 
// M/i + with j /qxj^3 Qs^/Jj'dxdv = is another solution such that, for /? > 4 

\\{v)^9s\\oo + \{vfgs\oo<^l, 

then fs = gs- Furthermore, if 9{x) is continuous on then Fg is continuous away from 
D. In particular, if Q is convex then T) = jq. On the other hand, if is not convex, 
we can construct a continuous function 'd{x) on in il.25\) with l^^loo ^ 1 such that the 
corresponding solution Fg in not continuous. 



We stress that the solution Fg is a genuine non equilibrium steady solution. Indeed, it is 
not a local Maxwellian because it would not satisfy equation (jO.ip . nor a global Maxwellian 
because it would not satisfy the boundary condition ()0.2p . 

We also remark that in addition to the wellposedness property, the continuity property in 
this theorem is the first step to understand higher regularity of Fg in a convex domain. The 
robust estimates used in the proof enable us to establish the following 5— expansion of 
Fg, which is crucial for deriving the necessary condition of the Fourier law (jl.ip . In the rest 
of this paper we assume the normalization M = 1. 

Theorem 1.2 Let 

fis = fJ- + Sfii + 6'^fj.2 -\ , (1-18) 

with J fJ.id'j = for all i from (T7)), and set 



Fg = n + ^fg. 

Then there exist /i,/2, ...,/m-i with || (f )'^e^l''l^/j||oo ^ 1, for ^ ^ C, < \,f3 > A, such that the 
following 5 -expansion is valid 

fg = 5fi + 6^f2 + --- + 6"'fi, 
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with \\{v)^ei\-\'fi\\^<lforOSC< 

jq^, /? > 4. In particular, fi satisfi^ 



les 



(1.19) 



1I7- 



h{x,v')^/j4^){n{x) • v'}dv' + 

n(x)-v'>0 V A* 



We have the following dynamical stability result: 

Theorem 1.3 For every fixed ^ C < 4q^> > 4 , there exist A > and eq > 0, depending 
on 60, such that if 



{vfe<\^\'[fiO) - fMoo + \{vfe<\^\'[f{0) - /,]U < eo, 



(1.20) 



then there exists a unique non-negative solution F{t) = fi + ^/Jlf{t) > to the dynamical 
problem and lll.5\) such that 



{vfe<\^\'[f{t)-moo + \{vfe'^\^\'[f{t)-f. 



< 



-xt 



If the domain is convex, 9{x) is continuous on dQ and moreover Fq{x,v) is continuous away 
from 7o and satisfies the compatibility condition 



Fo{x,v) = fi%x,v) 



n{x)-v' >0 



FQ{x,v'){n{x) ■ v')dv', 



then F{t,x,v) is continuous away from 79. 



The asymptotic stability of Fg further justifies the physical importance of such a steady 
state solution. We remember that, when 6 > 0, then Fg fails to be a global Maxwellian or 
even a local Maxwellian, and its stability analysis marks a drastic departure from relative 
entropy approach (e.g. [H]). Moreover, such an asymptotic stability plays a crucial role in 
our proof of non- negativity of Fg . 

An important consequence of Theorem 11.21 is the Corollary below which specializes the 
result to the case of a slab — ^ < x < ^ between two parallel plates kept at temperatures 
0(-i) = l-5and0(i) = l + (5. 

Corollary 1.4 Let x ^ ft = [—^, ^] and fs be the solution to 



vAfs + -^Lfs = -^T{fs, fs), 



1 1 

2<^<2' 



1 



yiJiv)fs{-,v) 



2tt[1 - 5]2 
1 



exp 
exp 



2[1 - 6] 

|2 



Ml<0 



\v\ 



2[1 + 6] 



ui^ li{u)fs{\,u)du, 



(1.21) 



uiy^ fi{u)fsi--,u)du, vi > 0, 



vi < 0. 



ui>0 



where vi is the component in the direction x of the velocity v. 
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Then fi, the first order in 6 correction to jj, according to the expansion of Theorem 
is the unique solution to 

vidJi + ^Lh = 0, -i<x<^, (1.22) 

1 i\ M ( ^ ^) 



^2' 



/i(o)^') = VM^^ / hii^.u)^/ ii{u){ui}du H 7==, < 



«i>o 2 •y//i(f 



In order to establish a criterion of validity of the Fourier law, let us remember that the 
temperature associated to the stationary solution Fg is given by 

Os{x) = ^ [ \v-Us\'^Fs{x,v)dv, (1.23) 
-aPs Jr3 

where ps = J-^a Fx{x,v)dv and Us{x) = pj^ f-^^ vFs{x,v)dv. The heat flux associated to the 
distribution is the vector field defined as 



Qs{x) = \i {v -Us(x))\v -Us{x)\^Fs{x,v)dv. (1.24) 
/r3 
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We state the Fourier law in the following formulation: there is a positive function k{6), 
the heat conductivity, such that U.l\) is valid for Fs . 

Theorem 1.5 // the Fourier law holds for Fs, then Oi{x) = 5 /j^s \v — uip/i is a linear 
function over —^<x<^. 

Available numerical simulations, Figure 1, [29] indicate the linearity is clearly violated for 
all finite Knudsen number Kn (= k in Figure 1). 

We therefore predict that the general Fourier law (jl.ip is invalid and inadequate in the 
kinetic regime with finite Knudsen number Kn. It is necessary to at least solve the linear 
Boltzmann equation (jl.l9p to capture the heat transfer in the Boltzmann theory. 

Without loss of generality, we may assume ^0 = 1 and Kn = 1 throughout the rest of the 
paper. Therefore in the rest of this paper we assume 

e{x) = l + 6^{x). (1.25) 

The key difficulty in the study of the steady Boltzmann equation lies in the fact that the 
usual entropic estimate for J FlnF, coming from dtF, is absent. The only a-priori estimate is 
given by the entropy production J Q{F, F) In F, which is very hard to use |4]. In the context of 
small perturbation of p, only the linearized dissipation rate || (I— P)/||^ is controlled0, and the 
key is to estimate the missing hydrodynamic part P/, in term of (I — P)/. This is a well-known 
basic question in the Boltzmann theory. Motivated by the studies of collisions in a plasma, 
a new nonlinear energy method in high Sobolev norms was initiated in the Boltzmann study, 
to estimate P/ ([H]) in terms of (I — P)/. For F = p + ^J^f , remembering the definitions 



^We denote by P/ the projection of / on the null space of L. 
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Figure 1: Nonlinearity of 9i in a kinetic regime 



(I1.14P of the linearized Boltzmann operator and (jl.l5p of the quadratic nonhnear collision 
operator, the Boltzmann equation (jl.2p can be rewritten for the perturbation / as 

{dt + v-V^ + L}f = r{f,f), (1.26) 

// the operator L were positive definite, then global solutions (for small /) could be con- 
structed "easily" for (jl.26p . However, L is only semi-positive 

(L/,/)>||{I-P}/||2, (1.27) 
where || • is the z^-weighted norm. The kernel (the hydrodynamic part) is given by 

P/ = {af{t, x)+v bf{t, x) + M_zi) c/(t, x)}^. (1.28) 

Note that we use slight different definitions of af{t,x) and Cf{t,x) from |19) to capture the 
crucial total mass constraint. The so-called 'hydrodynamic part' of /, P/, is the L^-projection 
on the kernel of L, for every given x. The novelty of such energy method is to show that L is 
indeed positive definite for small solutions / to the nonlinear Boltzmann equation (ll.26p . The 
key macroscopic equations for = (aj, bf, cj) connect P/ and {I — P}/ via the Boltzmann 
equation as in page 621 of [19] : 

A^/ ~ - P}/ + higher order terms, (1.29) 

where denotes some second order differential operator in [19]. Such hidden ellipticity in 
(k > 1) implies that the hydrodynamic part P/, missing from the lower bound of L, can 
be controlled via the microscopic part {I — P}/, so that L can control the full ||/||^ from 
(|1.27p . Such nonlinear energy framework has led to resolutions to several open problems in 
the Boltzmann theory [T9 | [20 | [25] . 
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It should be noted that ah of these results deal with idealized periodic domains, in which 
the solutions can remain smooth in for A; > 1. Of course, a gas is usually confined within 
a container, and its interaction with the boundary plays a crucial role both from physical 
and mathematical view points. Mathematically speaking, the phase boundary dVl x R'^ is 
always characteristic but not uniformly characteristic at the grazing set 70 = {(x, f ) : x G 
50, and v ■ n{x) = 0}. In particular, many of the natural physical boundary conditions 
create singularities in general domains ([22| I26j), for which the high Sobolev estimates break 
down in the crucial elliptic estimates (jl.29p . Discontinuities are expected to be created at the 
boundary, and then propagate inside a non-convex domain. Therefore completely new tools 
need to be developed. 

A new L'^-L°° framework is developed in \22\ (see [21j for a short summary of the method) 
to resolve such a difficulty in the Boltzmann theory, which leads to resolution of asymptotic 
stability of Maxwellians for specular reflection in an analytic and convex domain, and for 
diffuse reflection (with uniform temperature!) in general domains (no convexity is needed). 
We remark that the non-convex domains occur naturally for non-isothermal boundary (e.g. 
two non flat separated boundaries). Furthermore, the solutions to the boundary problems 
are shown to be continuous if the domain 17 is strictly convex. Different — L°° methods 
have been used in [21 |3l (SJ [I6l [17] in particular geometries. 

The new LP' — L°° framework introduced in [2^ has two parts: 

Positivity: Assume the wall temperature 6 is constant in (jl.Sp . It suffices to establish 
the following finite-time estimate 

\\Pg{s)\\lds<^j^ ||{I - P}5(s) 11^ boundary contributions! . (1.30) 

The natural attempt is to establish estimate for from the macroscopic equation (jl.29p . 
However, this is very challenging due to the fact that only / has trace in the sense of Green's 
identity (Lemma 12. 2p . neither P/ nor {I — P}/ even make sense on the boundary of a 
general domain. Instead, the proof of (jl.30p given in [22] is based on a delicate contradiction 
argument because it is difficult to estimate P/ via {I — P}/ directly in a setting in the 
elliptic equation (jl.29p . in the presence of boundary conditions. The heart of this argument, 
lies in an exact computation of P/ which leads to the contradiction. As a result, such an 
indirect method fails to provide a constructive estimate of (jl.SOp with explicit constants. 

L°° Bound: The method of characteristics can bootstrap the bound into a point-wise 
bound to close the nonlinear estimate. Let U{t) be the semigroup generated by v ■ V x + L 
and G{t) the semigroup generated by v ■ -|- i^, with the prescribed boundary conditions. 
By two iterations, one can establish: 

U{t) = G{t) + / G{t- si)KG{si)dsi + / / G{t - si)KG{si - s)KU{s)dsdsi. (1.31) 
Jo Jo Jo 

From the compactness property of K, the main contribution in (ll.3ip is roughly 

/ / \f{s,Xci{s;si,Xci{si;t,x,v),v'),v")\dv'dv"dsdsi. (1.32) 

Jo t)"bounded 

where Xci{s;t, x,v) denotes the generalized characteristics associated with specific boundary 
condition. A change of variable from v' to ^cl(") would transform the v' and ^''-integration 
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in (|1.32p into x and v integral of /, which decays from the theory. The key is to check if 



rXei(.;.i,Xei(.i;t,rr,.),.0|_^^^ (1.33) 
[ dv' J 

is vaUd. Without boundary, Xci(s; si, t, a;, w), w') is simply x — {t — si)v — (si — s)v' , 

and det | ^^^P | / most of the time. For specular or diffuse reflections, each type of 

characteristic trajectories Xc\ repeatedly interact with the boundary. To justify (|1.33p . various 
delicate arguments were invented to overcome different difficulties, and analytic boundary and 
convexity are needed for the specular case. 

Since its inception, this new — L°° approach has already led to new results in the study 
of relativistic Boltzmann equation [33], in hydrodynamic limits of the Boltzmann theory ([23] 
[23] [23]), in stability in the presence of a large external field ([15] [27]). 

Our current study of non-isothermal boundary {6 is non-uniform) is naturally based on 
such a — L°° framework. The main new difficulty in contrast to [22], however, is that the 
presence of a non-constant temperature creates non-homogeneous terms in both the linear 
Boltzmann (steady and unsteady) as well as in the boundary condition, so that the exact 
computation, crucial to the estimate (jl.30p . breaks down, and the L? — L°° scheme |22j 
collapses. 

The main technical advance in this paper is the development of a direct (constructive) 
and robust approach to establish (jl.30p in the presence of a non-uniform temperature diffuse 
boundary condition (j0.2p . Instead of using these macrosopic equations (e.g. ()1.29p ). whose 
own meaning is doubtful in a bounded domain, we resort to the basic Green's identity for the 
transport equation and choose proper test functions to recover ellipticity estimates for o, 6, c 
and hence P/ directly. In light of the energy identity, {1 — -Py}/ is controlled at the boundary 
7+, but not Py/ in ()1.16p . The essence of the method is to choose a test function which can 
eliminate the P^f contribution at the boundary, and to control the a, 6, c component of P/ 
respectively in the bulk at the same time. The choice of the test function for c is rather direct: 
we set 

for some constant (3c to be determined, with —A(f)c = c and (pc = on dfl. On the other 
hand, the test function for b is rather delicate. In fact, two different sets of test functions 
have to be constructed: 

{vf - ^b)VJj'dj(l)l, i,j = l,...d, 
\v\'^ViVj,/JIdjcl)l, i / j, 

with —^(pl = bj and = on d^l. In particular, a unique constant /3c can be chosen to 
deduce the estimates for b, thanks to the special structure of the transport equation and the 
diffuse boundary condition. The choice of the test function for a requires special attention. It 
turns out that in order to eliminate the contribution P^/, we need to choose the test function 

(|i;|^ - (3a)v ■ Vx(pa^/Jj; 

with — A(/>a = a and Neumann boundary condition 

d 

—(pa = on dn. 
on 
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This is only possible if the total mass of / is zero, i.e., 

// f{x,v)./JIdxdv = V2^ a{x)dx = 0. (1.34) 

This illustrates the importance of the mass constraint, which unfortunately is not valid for 
the steady problem (jO.ip and (j0.2p . So we are forced to use a penalization procedure (see 
below) to deal with it. The key lemma which delivers the basic estimates is Lemma 13.41 
Furthermore, in the dynamical case, also the time derivatives dt4>c, dt(j)l and dtcpa need to 
be controlled in negative Sobolev spaces. This is possible due to the special structure of the 
Boltzmann equation as well as the diffuse boundary condition. Once again, the total mass 
zero condition (|1.34p is essential. This is the key to prove the crucial Lemma [6. 2 1 Even though 
this new unified procedure can be viewed as a 'weak version' of the macroscopic equations 
(|1.29p . the estimates we obtain via this approach are more general. For instance, the p.29p 
was only valid for dimension > 2, but the new estimates are valid for any dimension. There 
seems to be a very rich structure in the linear Boltzmann equation. 

To bootstrap such a estimate into a L°° estimate, we define the stochastic cycles for 
the generalized characteristic lines interacting with the boundary: 

Definition 1.6 (Stochastic Cycles) Fixed any point {t,x,v) with{x,v) ^70, let {tQ,XQ,V()) 
= {t,x,v). For Vk+i such that Vk+i-n{xk+i) > 0, define the (k + l) -component of the back-time 
cycle as 

{tk+i,Xk+i,Vk+i) = (tk - th{xk,Vk),Xb{xk,Vk),Vk+i). (1.35) 

Set 

Xc\{s;t,x,v) = ^ l[t^^^_t^)(s){xfc + (s - tfc)!;fc}, 

k 

Vci{s;t,x,v) = J2Mtk+i,tk)(^)'"k, Vci{s;t,x,v) = {Vci{s;t,x,v),v). 

k 

Define Vk+i = {v £ | v ■ n{xk+i) > 0}, and let the iterated integral for k > 2 be defined as 





[1.. 







dai, (1.36) 



where daj = fi{v){n{xj) ■ v)dv is a probability measure. 



We note that the Uj's (j = 1,2,...) are all independent variables, and tk,Xk depend on 
tj,Xj,Vj for j < k—1. However, the phase space Vj implicitly depends on (t, x, v, vi,V2, ...Vj-i). 
Our method is to use the Vidav's two iterations argument (^37J) and estimate the L°°-norm 
along these stochastic cycles with corresponding phase spaces H^zlVj. The key is to esti- 
mate measures of various sets in n'^z|V,- in Lemma l4.2i We designed an abstract and unified 
iteration (14. Sp . which is suitable for both steady and unsteady cases. New precise estimate 
(|4.16p of non-homogeneous terms resulting from the non-constant temperature are obtained 
in Prop. 14.11 Based on a delicate change of variables in Lemma 12. 3^ such an estimate is 
crucial in the proof of formation of singularity for a non-convex domain in Theorem 11.11 

In order to keep the mass zero condition and to start iterating scheme, it is essential to 
introduce a penalization e to solve the problem 

e/+' + t;-V,/+i + L/+i = r(/,/). 
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The presence of e ensures the critical zero mass condition (jl.34p . It is important to note that 
our L°° estimate are intertwined with at every step of the approximations, which ensures 
the preservation of the continuity for a convex domain. The continuity properties of our final 
solutions follows from the limit at every step. Moreover, the proof of continuity away 
from the singular set T) in (jl.lip in a general domain is consequence of a delicate result for 
Qgain in 126]. 

To illustrate the subtle nature of our construction, we remark that for the natural positivi- 
ty-preserving scheme: 



we are unable to prove the convergence, due to breakdown of L/ , whence the mass zero 
constraint ()1.34p fails to be satisfied. Consequently, we are unable to prove > in our 
construction. Such a positivity is only proven via the dynamical asymptotic stability of -F^, 
in which the initial positivity plus the choice of a small time interval are crucial to guarantee 
the convergence of the analog of (I1.37P in the dynamical setting. 

Our estimates are robust and allow us to expand our steady state Fs in terms of 5, the 
magnitude of the perturbation. This leads to the first order precise characterization of Fg by 
/i, which satisfies (I1.19p . 

It should be pointed out the our rather complete study of the non-isothermal boundary 
for the Boltzmann theory for 5 <C 1 forms a mathematically solid foundation for influential 
work in applied physics and engineering such as in [31)|32]. where the existence of such 
steady solutions to the Boltzmann equation is a starting point, but without mathematical 
justification. We expect our solutions as well as our new estimates would lead to many new 
developments along this direction. 

The plan of the paper is the following. In next section we present some background 
material and in particular a version of the Ukai Trace Theorem and the Green identity as well 
as a new estimate at the boundary in Lemma 2.3. Section 3 is devoted to the construction 
of solutions to the stationary linearized problem. In particular, we prove Lemma [3.41 which 
provides the basic estimate of the norms of P/ in terms of the norm of (I — P)/. In 
Section 4, after introducing an abstract iteration scheme suitable for proving L°° bounds, we 
prove, in Proposition 14. !( the existence of the solution to the linearized problem in L°° and 
related bounds. In Section 5 we combine the results of the previous sections to construct 
the stationary solution to the Boltzmann equation and discuss its regularity properties. In 
particular, we give the proof of the 5-expansion and use it to establish a necessary condition for 
the validity of the Fourier law. In Section 6 we extend the estimates to the time dependent 
problem. Section 7 contains the extension of the L°° estimates to the time dependent problem 
and the proof of the exponential asymptotic stability of the stationary solution. From this 
we then obtain its positivity. 

2 Background 

In this section we state basic preliminaries. First we shall clarify the notations of functional 
spaces and norms: we use || • ||p for both of the L^{Q x R^) norm and the L^{Q) norm, and ( • , • ) 




■l+l 



(1.37) 
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for the standard L^(r2 X R^) inner product. Moreover we denote = ||i^^''^-||2 and = 
II/II2 + X^i=i ll^x/lb- For the phase boundary integration, we define d'y = \n{x) ■ v\dS{x)dv 
where dS{x) is the surface measure and define |/|p = \ f{x,v)Y'd'y = \f{x,v)\P and the 
corresponding space as LP{dQ. x R^;(i7) = L'P{d^ x R^). Further |/|p,± = |/l-y^|p and 
l/|oo,± = |/l7±U- We also use \f\l = /^^ |/(x)| = !Q^\f{x)\P. Denote /± = 

Recah that x G C R'^ for d = 1,2, 3, and v G R'^. For A S R'^, the notation v ■ A means 

In the next lemma, we prove that the trace of / is well-defined locally for a certain class 
of/: 

Lemma 2.1 (Ukai Trace Theorem) Define 

7^ = {{x,v) G 7 : \n{x) ■ v\ >e, \v\ < -}. (2.1) 

Then 

l/Vli <e,n ll/l|i + ll^-V,}/||i, (steady) 

|/V(r)|idr S,n / \\f{r)\\idT+f + ^ • Vj/(r)||idr, 
J s J s 

(dynamic) 

for aUO<s<t. 

Proof. Recall the notation (jl.8l) and the definition of backward exit time t]j(x,v) = t\y(x,v) 
and Xb(x, v) = Xh(x, v) = x — t\3(x, v)v S dQ in (jl.9p . 

In the steady case, from [10], page 247, we have the following identity: 



f(x,v)dxdv = / / f(x — sv,v)dsd'j 



+ / f(x + sv,v)dsdj. 



'7- JO 

Assume ||/||i, \\{v • Vx}f\\i < 00. Then f(x + su,f) is an absolutely continuous function of 
s for fixed (x,v), so that, by the fundamental theorem of calculus 

f(x,v) = f(x - s'v,v) + J V ■Vxf(x + TV,v)dT 

for (s ,x,v) e [0,tb(a;,f)] x 7+, 

r' 

f(x,v) = f(x + s'v,v)+ / f • Va;/(x + Tf , u)(ir 

JO 

for (s',x,t>) E [0,tb(a;,-^^)] X 7_, 

and therefore, for both cases, 



\f(x,v)\<\f(x-s'v,v)\+ \vV,f(x-Tv,v)\dr, (2.2) 

\f(x, V)\ < \f(x + S'V, V)\+ [''^ ' V • ^xf(x + TV, v)\dT. (2.3) 
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On the other hand, for x S dVl assume that is locally parameterized by : {y G R*^ : 
|rE — ?/| < (5} — R so that 

\{x — y) ■ nix)\ ,„9 ^, 

sup I " ^' < max \Vli{y)\. 

\x-y\<5 \x-y\<S 

By the compactness of (and Sfi), uniformly in x, we have \{x — y) ■ n{x)\ < Cfi\x — for 
all y G dQ. Taking inner product of x — Xb(x, v) = t^{x, v)v with n{x), we get 

th{x,v)\v ■ n{x)\ = \{x- Xhix,v)) ■ n{x)\ < Cn\x - Xb(x,f)|^ = Cn\v\'^\tb{x,v)\'^ . 

Therefore we deduce 



tb{x,v) > ■ (2.4) 



\n{x) ■ V 

If {x, v) G (1 7+ or (x, v) G n 7_ then, by the definition of 7*^ and (40) in 

mm{ti,(x,v), tb(x,-w)} > ^"^'^i "i!"* > Cf^V. 

First we integrate (|2.2p . for (x,f) G 7^07+, and then for s' £ [0,tY,{x,v)]. Similarly integrate 
(|2.3p . for {x,v) G 7^ n 7_, and then for s' G [0,tb(a^5 — ^)] to have 

C^V|/V|i<C{||/||i + ||t;-V,/||i}. 

In the dynamic case, from |1U] in page 247, we have the following identity : 

f{s,x,v)dsdxdv = / / f{s,x—{t — s)v,v)dsd'j 

[0,t]xnxR3 J-y+ Jmax{0,t-tb{x,v)} 

p i-mm{t,ti,ix,~v)} 

+ / f{s,x + sv,v)dsd'y. (2.5) 



J7_ Jo 

Assume again ||/||i, + v • Va;}/||i < 00. Tthen f{s,x + sv,v) is absolutely continuous of 
s for fixed x,v, so that, by the fundamental theorem of calculus 

f{s,x,v) = f{s-s',x-s'v,v)+ vVxf{T,x - {s - T)v,v)dT 

Js-s' 

for (s, s ,x, v) G [0, t] X [max{0, s — t\^{x, v)}, s\ x 7+, 

/s+s' 
V ■ ^xf{T, X + TV, v)dT 

for (s, s', X, v) G [0, t\ X [0, min{t — s, tb(3^) ~^)}] ^ 7-; 

and hence, for both cases, 

|/(s,x,u)| < |/(s - s',x - s'z;,u)| + / I?; • V^./(r,x - (s - T)u,T;)|dr, 



J s—s' 

/s+s' 
\V ■ VxfiT,X + TV,v)\dT. 

Then the rest part of proof is exactly the same as in the steady case. ■ 
The following Green's identities are important in this paper. 
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Lemma 2.2 For the steady case, assume that f{x,v), g{x,v) G L^(ilxR^), v-Vxf, f •V^.g G 
L'^{n X R3) and /^,c/^ G ^2(017 x R^). T/ien 

ff {v VJ}g + {v ■ V^g}/ dvdx = [ fgdj. (2.6) 

For the dynamic case, assume that f{t,x,v), g{t,x,v) £ L'^{[0,T]; L'^{Q,x'R'^)), dtf + v-'Vxf , 
dtg + v ■ V^g G L^{[0,T] x n x K^) and f-y,gy G L^{[0,T] x dQ x R^). Then, for almost all 

t,se[o,T], 

t ff {dtf + V ■ Vxf}9 dvdx + f' ff {dtg + v- V^g}! dvdx 

Js JjQxR^ Js JJnx'R^ 

= ff f{t)g{t)dvdx - ff f{s)g{s)dvdx + f f fgd^dr. (2.7) 

JJnxR'i JJnxR^ Js J-y 

Proof. See the proof in Chapter 9 of [lOj. ■ 

Lemma 2.3 Let C R'^ and recall the notation il.8\) . If 

V ■ n(xb(a;, v)) < 0, 
then t\3{x,v) and Xi,{x,v) are smooth functions of {x,v) so that 

n{xi,) ti,n{x\3) 

Vx^b = 7 — r> Vtjtb = —- 7 — r> (2.8) 

V ■ n[xb) V ■ n(xbj 

VxXh = I - V^tb O V, VyXh = -tbi - Vi;tb ® V. 
Ford = 2,3, if x e dQ 

f \f{x-tb{x,u)u\\n{x)-u\du< f \f{y)\dS{y), (2.9) 
ys'*-i Jdn 

and ifxGQ and {y G R'' : \x - y\ < e|} n = 



|/(x-tb(x,u)u)|duS / \f{y)\dS{y), (2.10) 
S'*-! Jan 

Proof. Assume that dft is locally parameterized by ^ : R'^ — )• R and Q is locally {x G R'^ : 
^ (x) < 0} so that 

^{xb{t,x,v)) = C{x - tb{x,v)v) = 0. 
By the implicit function theorem, taking Xj and vj derivative respectively, 

Oji - 2^ Oki—Vk = 0, so that 



k-i ' ^■''^(^^b)' 

dtb_ „ , dtb , nj{xb) 



({jCtb + > , OkK^Vk = 0, so that — = -tb^ . ^ ■ 
ovj ovj V ■ n(Xb) 



k=l 

We then have 



'9(xb)i _ ^ t'ira-,(xb) tj(xb)i _ Vin-jixb) \ 
dxj V ■ n{xb) ' dvj \ v ■ n{xb) J 
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Now we prove (j2.9p for d = 3. Without loss of generality we may assume that 93^(2;b) 7^ 0. 
Using the spherical coordinates, u = (sin 6 cos sin 9 sin cos 6) £ and 

^(xi — tb(2;, u) sin 6 cos 2:2 — tb(2;, u) sin cos </>, 2:3 — tb(2;) u) cos 0) = 0. 

By the implicit function theorem we compute 

sin 9 + 9i^(xb) cos 9 cos (/> + d2C{xh) cos sine/) 



— thix,u{9,(j))) = -tb{x,u{9,(l))) MV7 M 

5(9 u(6',</>) •n(xb)|V^C(2;b)| 

d ro w ,^s9i^(xb)sin6'sin(/) - (92C(xb)sin6'cos0 



dcP"' ' ' ' ^ u(e,</.)-n(xb)|V,e(xb)| 

Since the Jacobian matrix is 



Jac 



[ g((a^b)i,(xb)2) 



we have 



and hence 



- sin 9 cos — tb cos cos (j) — sin 6 cos + tb sin sin { 

- sin 9 sin (/> — tb cos sin (j) — sin 9 sin — tb sin cos ( 



/ a((xb)i,(xb)2) \ (tb)^sing|a3g(xb)| , . 

V 5(0,0) y - |n(xb)-u||VC(xb)|' ^ ■ ^ 



|/(x-tb(x,u')u')||n(x)-u'|du' 

^ / /'G[o.) </,'e[0 2.) - *b(2;,u')u')||n(x) • u'| sin0'#'d0', (2.12) 

(e'',0')~{0,<^) 

where u' = (sin 9' cos 0', sin 9' sin </)', cos 9'). Now we apply a change of variables 

{9',(p') ^ ((xb(x,u'))i,(xb(x,u'))2) ^ ((x'b)i,(x'b)2), (2.13) 
and use ()2.1ip to further bound (|2.12p as 

|tb(x,u(x;))p |a3C(x'b)|'^^"^^^^^"^^^- ^'-''^ 

(X(^)2~(Xb)2 

|V^(x' )| 

Notice the surface measure of dVt is dS = \^) Y^ d{x'^)id{x'^)2 if d-ii{x'^) / 0. Denote 
^ Ic'sUa^b)! 

11(2^11) = r~ — "TT ^-iid we use (12. 4D as 
|x — Xjjl 

tb(2;,u(xb)) >n \n{x) •u(xb)|, tb(2;,u(xb)) >n |n(xb) • u(xb)|, 
to bound (j2H]) by 

\f{y)\dS{y). 
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For d = 2, using u = (cos 6*, sin 6*) if c?2^(2;b) 7^ we can compute 
d , , n(xb) • (sin0, — cos^) 



dd ' ' ' ' " n(xb) -^0) 

and 

, , ( dMx \ ^ tb(x,u(g))|92C(xb)| 

"""'V 50 ) |n(xb)-u(0)||Ve(xb)r ^ 

The rest of proof of (j2.9p is same (even simpler!) as the d = 3 case. 

For (j2.10p . since there is a lower bound for t-^ix^'ti) > e for {\x — y\ < e} n 50 = and 
u e S°'~\ using pTTjl and (l2l^ . it is easy to prove (j2TT(l . ■ 



3 Estimate 

The main purpose of this section is to prove the following: 
Proposition 3.1 Assume 

il g{x,v)^/Ji dxdv = I r^yjld'y = 0. (3.1) 

J JflxR^ J-y- 
Then there exists a unique solution to 

vV,f + Lf = g, /_ = P^f + r, (3.2) 
such that //f^xR3 fV^ dxdv = and 

+ I/I2 < II5II2 + \r\2- 



For the proof of Proposition 13 . 1 1 we need several lemmas. We start with the simple transport 
equation with a penalization term: 

Lemma 3.2 For any e > 0, there exists a unique solution to 

ef + v- V^./ = g, /_ = r, 

so that 



12 + I/I2 S 115112 + kb, 

for all (3 > Q^C, > 0. Moreover if g and r are continuous away from the grazing set 70, then 
f is continuous away from D . In particular, if Q is convex then 2) = 70 . 
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Proof. The existence of / and L'^-bound follow from integration along the characteristic 
lines of ^ = v e R"', and ^ = (Recall (jl.Sp for notations). More precisely, setting 
h{x,v) = {v)^e''^^^' f{x,v), the integrated form of the equation for h is: 



h{x,v) = ltyt^(^^^y)^h{x -tv,v)e ^* 

t 



+li<ib{x,^) { (^^>''e'^l^l V(x - tv, v)e 
ft 



-et 



+ 



{v)^e'^^''^'g{x -{t- s)v,v)e-<'-'Usy 



where t\3{x,v) is defined in (|1.9p . We prove the L°°— bound by choosing a large t = t{e) such 
that 

\h{x,v)\ < ^{\\h\\^ + \h\^} + \\{v)^e<\^\\\\^ + e-'\\{vfe<\^\'g\U 

In order to prove the continuity, let {x,v) £ Q x R^\X'. Then by the definition of D, 
n(xb(x,u)) ■ V < and hence t\3(x,v) is smooth by Lemma 12.31 Therefore, if g and r are 
continuous, f{x,v) is continuous at {x,v) £ fl x 'R'^\D. 

Suppose now that Q, is convex. In order to show that D = jq, since D 70 is true 
for any 17, it suffices to show that 2) C 70- Since Q is convex t\3{x,v) = = tiy{x,—v) for 
X S d^l, V / 0. Therefore 79 = dft x {0}. Hence, if {xb{x,v),v) G 7q then v = and 
X]3{x,0) = X G dQ and {x,v) G 7q . 

The L^-estimate and the uniqueness follow from Green's formula, since ||/||2 < ||(t')''e^l'^l /||oo- 



In next lemma we add to the penalized transport equation a suitably cut-offed linearized 
Boltzmann operator. Moreover we include a reduced diffuse reflection boundary condition, 
with the purpose of setting up a contracting map argument. We have 

Lemma 3.3 For any e > 0, m > 0, and for any integer j > 0, there exists a unique solution 
to 

ef + v-VJ + L^f = g, /- = (!- -)P^f + r, (3.3) 

J 

with Lm the linearized Boltzmann operator corresponding to the cut-offed cross section Bm = 
min{i?, m}. Moreover, uniformly in j, we have 



\u + I/I2 <e,m hh + kb- 

Finally the limit f as j ^ 00 of the sequence {/•'} exists and solves uniquely 

ef + v-VJ + Lmf = g, f-= P^f + r. (3.4) 

Proof. Denote Lm = Vm — Km- For any j, we apply Lemma [3^2] to the following double- 
iteration in both j and i: 

ef'^^ + v-V^f'+^ + Vmf'^^-Kmf' = g, (3.5) 

J 
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with /° = 0. 

Step 1. Wc first fix m and j and take ^ — >■ oo. 
From Green's identity, 

^ J 

Prom {Kmif + f+'),f + f+^) < {vmif + + /^+'), we deduce 

Moreover, there is Cj such that 

1(1 - \)P,f + r\l_ < 1(1 - \)P,f\l. + ^\P,f\l- + CMl (3.6) 
3 J 



We derive from \P-s,j%,- < and || • lb > ^ 



l\2 ^ I fi\2 II . ||„ > J- 

11 l'-^ — TO 

1 1 „tf I 1 lO - 1 11 „tf iiO 1 1 



f^+l||2 



Since ^ + ^ — ^ > ^ and 1 — j + < 1, by iteration over £, for some r/g^^j- < 1, 

Wf^'WL + < VcmAWfWL + + Ce,mA\r\l + llfflll}- 

Taking the difference of /^+^ — we conclude that is a Cauchy sequence. We take i ^ 00 
to obtain as a solution to the equation 

ef + • V./^' + Lmf = 9, f- = {l- j)Pjf' + r. (3.7) 

Step 2. We take j ^ 00 for 

By Green's identity we obtain uniformly in 

ewpwl + (w^/^■) + ^1/^11+ - ^1^7/^+^11- = / /V 

We rewrite for any 77 > 

\\P,P+r\l- = \\P^P\l- + \\r\l + J^P-yfrd^ 

< l\Pyf\l- + Crj\r\l + v\P^f\l-, (3.8) 
so that from J f^g < |||/^||2 + C'slls'lli and from the spectral gap of L^, we have 



WfWl + 11(1 - P)f\\l^ + i|(l - P,)f\l^ < C,A\r\l + Ml} + v\P,f\l- (3.9) 
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But from the equation we have 

V ■ Vx.[/^f = -2e[ff - 2fLm{I - F)f + 2fg. 

Taking absolute value and integrating on x R'^, from (j3.9p we have 

||^;-V.(/^yi|i < Cs{\\f\\l + Ul-P)P\\t + \\9\\l} 
< C^A\r\l + \\9f2} + vCe\P-yf\l-. 

Hence, by Lemma l2.1i for any 7"^' in (|2.ip away from 70, we have 

\P Vli < C,,r,AMl + II5II2} + r]C,,,,\P^P\l_. (3.10) 

From (jl.lGp we can write P^f^ = z^{x)^^ for a suitable function z^{x) and, from \P^f^l^^i\2 < 
I/-' l-^e'b < +00, for e' small 



7' 

|2 



P^/-^ 1 ./I2 = / \z^{x)\'^ j fj,{v)\v ■ n{x)\dvdx 

Jan ■^|n(x)-u|>£'>|<^ 

> / \z.y{x)\'^dx X - 1 fi{v)\v ■ n{x)\dv 



\\P,P\l (3.11) 



where we used the fact 



e 2 



li{v)\n{x) ■ v\dv < e''"\\ \v\\\dv\\ e'l^^l /^dT;^ < Ce', (3.12) 

|n(z)-i;|<£' J-e' Ju.^ 

fi{v)\n{x) ■ v\dv < Ce' . 

\v\>l/e' 



Therefore we conclude 



\\p^P\l-\{^-p,)P\l+ < |p^/-''vii-l(i-^7)/'Vl2,+ (3-13) 



< C,^v,e'{\r\l + \\9\\l} + vCe,e'\Pjf\l-. 

Adding 4 x I^M) to (f3l^ . we obtain: 

25||/^||i + 1(1 - P,)/^||+ + 4||(I - P)f\\l^ + \\P,P\l 

< Ce,,{\r\l + + r?(l + C,,,OI^'7/'li,- 

Choosing r] small and taking the weak limit j — )• 00, we complete the proof of the lemma. ■ 

Next lemma states the crucial bound for P/. It will provide uniform in e estimates 
which allow to take the limit e — )• in (13.41). 
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Lemma 3.4 Let f be a solution, in the sense of i3.15\} below, to 

vVJ + Lf = g, f_ = P^f + r, (3.14) 

with 

/ / fy/Jidxdv = gy/Jidxdv = / r^/JIdj = 0, 

then we have 

W^fWl < 11(1 - P)f\\l + Ml + 1(1 - P,)f\l+ + \r\l 
Proof. The Green's identity ()2.6p provides the following weak version of (j3.14p : 

[i^fdj-ff vV,ijf = -[[ vMi-P)/+// i^a- (3.15) 

J7 JJnxR^ JJnxR^ JJnxii^ 

Recall that P/ = {a+u •5+c[^^^ — flly^ on ^7 x R^. The key of the proof is to choose suitable 
test functions to estimate a, b and c in p.22p . (|3.29p . ()3.34p . p.42p thus concluding the 
proof of Proposition 13. li 

Step 1. Estimate of c 

To estimate c, we first choose the test function 

i; = i;c={\vf -^c)VJIv-VMx), (3.16) 

where 

-A^cpcix) = c{x), (j)c\dn = 0, 
and /3c is a constant to be determined. From the standard elliptic estimate, we have 

W'PcWh^ < ||c||2. 

With the choice ()3.16p . the right hand side of (j3.15p is bounded by 

||c||2{||(I-P)/||2 + |bl|2}. (3.17) 

We have 

d 

V ■ Vxipc = ^ (bp - /3c)viy/JIvjdij4>cix), 
so that the left hand side of (|3.15p takes the form, for i = 1, • • • ,d, 

d 



I (n(x) • t>)(|7;p - /3c)V)"y]^i^j</'c/ 
Janx'R? ~[ 

-II - Pc)^{ y\ ViVkdi,4>c]f. (3.18) 



We decompose 



= P^/ + l^+(l-P^)/ + l^_r, on 7, (3.19) 

/ = {a + t;.6 + c[^-|}^ + (I-P)/, onJ]xR3. (320) 
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We shall choose /3c so that, for all i 

{\v\^-l3c)vffiiv)dv = 0. (3.21) 



1 _ Ml 

Since = 2 the desired value of /3c is /3c = 5. Because of the choice of /3c, there is no 

a contribution in the bulk and no P^f contribution at the boundary in (|3.18p . 

Therefore, substituting (|3.19p and (j3.20p into (13.18p . since the b terms and the off-diagonal 
c terms also vanish by oddness in v in the bulk, the left hand side of (j3.15p becomes 



d 

-^2 - Pc)vK^ - i^)fJ-iv)dv dii(j)c{x)c{x)dx 

~^^R3 2 2 Jq 

-Y, ff (I^P - Pc)v^VJ^{v • V,.)9,0c(I - P)/. 

~l J JQxR^ 

For /3c = 5, /i^3(bP - /3c)vf - = 107rV2^. Therefore, we obtain from (pTTTll 

- W7tV2^ [ A,M^)c{x) < ||c||2{|(l - P^)f\2,+ + 11(1 - P)/||2 + II5II2 + |r|2}, 
Jn 

where we have used the elliptic estimate and the trace estimate: 

|Vx(Ac|2 < \\(l>c\\m ^ l|c||2- 

Since —Ax(j)c = c, from (|3.16p we obtain 

Ml < {1(1 -P^)/|2,+ + 11(1 -P)/||2 + ||5l|2 + |r|2}||c||2, 

and hence 

llclli < 1(1 - P,)f\l+ + 11(1 - P)/||i + ll^lli + \r\l (3.22) 
Step 2. Estimate of b 

We shall establish the estimate of b by estimating [didjA~^bj)bi for all i,j = l,...,d, and 
{djdjA-^bi)bi for i 7^ j. 

We fix To estimate didjA^^bjbi we choose as test function in ()3.15p 

^ = = (vf - ^)^d,<l,l i, j = 1, . . . , d, (3.23) 

where f5b is a constant to be determined, and 

- A^(l)l{x) = bj{x), (l)l\dn = 0. (3.24) 

Prom the standard elliptic estimate 

WIh- < Wbh. 
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Hence the right hand side of (|3.15p is now bounded by 

||6||2{||(I-P)/||2 + ||5||2}. (3.25) 

Now substitute (IXTOll and IK20\\ into the left hand side of (IXTKI) . Note that {vf - (3b){n{x)-v}n 
is odd in v, therefore P^/ contributions to ()3.15p vanishes. Moreover, by ()3.20p . the a,c 
contributions to ()3.15p also vanish by oddness. Therefore the left hand side of ()3.15p takes 
the form 



(n(x) • v){vf - /3b)VJldAf -II i^f - fib)^{Y, vidijcPUf 
= I {n{x)-v){vf-l3b)Vlidj4[{l-P^)f + r]l^^ (3.26) 

-II yZivf - f3b)vffidij(l)l{x)bidvdx (3.27) 
Jn JR3 I 

~ II E(^«' - Mvi^di,^l{x){I - P)/. 

|2 

1 tt3 

Furthermore, since fx{v) = 27 ni=i ^ ~ '^^^ choose /3f, > such that for all i, 

I [{v^f - PbMv)dv = I [vl - Pb]e--^dvi = 0. (3.28) 

We remark that the choice (|3.28p also plays a crucial rule in the dynamical estimate (|6.1ip . 

Since ^{v) = i^e 2 , the desired value is = 1. 

Note that for such chosen Pi,, and for z 7^ fc, by an explicit computation 

9 9 f 9-1 l^ll'^ 1^21"^ |tJ3 1"^ 

- Pb)v^ndv = / (i^i - /3b)w2— e ~e 2 dv 



2-K 

9 it-ii^ 
{v{ - pb)e 2 dvi = 0, 



R 

{vf - (3b)vfndv = [vf- (3bvf]e-^dvi = 2^2^ / 0. 
Jn 

Therefore, (13:2711 becomes, by I^^M, 

- {vi-Mvffidvdij<f>l{x)bi+y2i i'^i - Mvlfi] dkj(t>i{x)bk 

JjnxR3 ,7T^ -^R-^ J^ 

k(^i) s ^ / 

=0 

= 2^2^ I {didjA-^bj)bi. 
Jn 

Hence we have the following estimate for all i,j, by (j3.25p : 

d,d,A-\b,^ < 1(1 - P,)f\l^ + 11(1 - P)/||i + ll^lli + \r\l+ e\\b\\l (3.29) 

In order to estimate dj{dj/Sr^bi)bi for i 7^ j, we choose as test function in (I3.15P 

= \v\^ViVj^dj4>l{x), i / j, (3.30) 
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where is given by (|3.24p . Clearly, the right hand side of (j3.15p is again bounded by (j3.25p . 
We substitute again (j3.19p and (|3.20p into the left hand side of (|3.15p . The Pjf contribution 
and o, c contributions vanish again due to oddness. Then the left hand side of (j3.15p becomes 



{n ■ v}\v\'^ViVj^/Jldj(plf - \v\'^ViVj^{S2vkdkj(l)l}f 



7+ 



\v\^vfvjid[dij(j)lbj + djj4>l{x)bi] 



-11 \v\'^ViVjVky/Jldkj4>l{x)\l-'P]f. 

Note that p.32p is evaluated as 

7V2^ [ {{didjA-\)bj + {djdjA-\)h}. 

Jn 

Furthermore, by \dj(t>l\2 < W\h'^ < Ph, so that 

(I33I1) + » < ||6||2{|(1 - P,)fk+ + \r\2 + 11(1 - P)/||2}. 
Combining (I3.29p . we have the following estimate for i ^ j, 

djdjlSr^bihi 



(3.31) 
(3.32) 
(3.33) 



< 



did,A-\bj 



+ 1(1 -^7)/li+ + 11(1 -P)/lli + 



\l + \r\l + e\ 



< 1(1 - P,)f\l+ + 11(1 - P)/lli + Ml + \r\l + e\\b\\l 
Moreover, by IK29\\ . for i = j = 1, 2, . . . , d, 

d^A-'bjbj 



(3.34) 



< 1(1 - P,)f\l+ + 11(1 - P)/ll^ + Ml + \r\l + e\\b\\l 



(3.35) 



Combining (I3.34|) and (I3.35p . we sum over j = 1, 2, . . . , d, to obtain, for alH = 1, 2, . . . , d. 



2< \ii-p,)f\i+ + m-F)f\\i + Ml + \r\i 

Step 3. Estimate of a 

The estimate for a is more delicate because it requires the zero mass condition 



(3.36) 



// 



adx. 



axR3 



fy/Jidxdv = = / 
Jo. 

We choose a test function 

d 



(3.37) 



i=l 
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where 

d 

-A^.(/)a(x) = a(x), ^</'a|9n = 0. 
It follows from the elliptic estimate with a = that we have 

W't'aWH-i < \\a\\2- 

Since /j^3(^^^^ lJi-{v)dv ^ 0, we can choose /3a > so that, for all i, 

I (H'-/5a)(^-|)K)VW=0. (3.38) 
7r3 2, Z 

Since ^{v) = 2~, the desired value is /3a = 10. Plugging ■(/;„ into p.lSp and its right 

hand side is again bounded by 

||a||2{||(I-P)/||2 + ||5ll2}. 

By ()3.19|) and (j3.20p . since the c contribution vanishes in (j3.15p due to our choice of /3a and 
the b contribution vanishes in (|3.15p due to the oddness, the right hand side of (|3.15p takes 
the form of 



^ / {n • v}{\v\^ - (3a)vi^diM^)[P^f + Pi)fh+ + ^^17+] (3-39) 
i=i "'t 

ff {\v\^ - Mv^Vk^^kMx)a{x)^^{v) (3-40) 

If {\V\'' - Pa)v^Vk^Mx){l-'P)f. (3.41) 



i,k=l 
d 

We make an orthogonal decomposition at the boundary, 

Vi = {v ■ n)ni + {v±)i = VnTii + {v±)i. 
The contribution of P^f = z.y{x)y/]I in (j3.39p is 

(|w|^ - /3a)v ■ Vx(pa{x)Vnf^iv)z^{x 
(bp - /3a)Vn^^VnHiv)z-y{x) 



7 



dn 

The crucial choice of (|3.38p makes the first term vanish due to the Neumann boundary con- 
dition, while the second term also vanishes due to the oddness of {v±)iVn for all i. Therefore, 
dOO]) and (tOT]) are bounded by 

||a||2{||(I-P)/||2 + |(l-P,)/|2,+ + |r|2}. 
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The second term (|3.40p . for k ^ i vanishes due to the oddness. Hence we only have the k = i 
contribution: 



1=1 

Using —/S.x<i>a = a we obtain 



yZ {\v\'^ - (3a)iVif fidii(j)aa. 



a 



\l < 11(1 - P)/lli + 1(1 - P-y)f\l+ + \r\l + Ml (3.42) 



We close this section by proving Proposition I3.1[ 

Proof of the Proposition 13.11 

We keep m fixed and take e — )• in Lemma [3.3l by using Lemma [3. 41 which obviously holds 
even with the additional penalization term. Indeed J /qxrs — ^ — ^\//^^^' H^nce we 
have e //qxrs f^\/T^ — therefore, for any e > 0, 



f^yjjldxdv = 0. 
We first obtain 

llPrili < 11(1 - P)rili + 1(1 - P,)r\l+ + \r\l + Ml + eWfWl (3.43) 
On the other hand from Green's identity: 

e\\f\\l + {Lmf,n + \\f\l+-\\P,f + r\l_= j fg, 

we deduce from the spectral gap of 

ewrwi + 11(1 - p)n\i^ + - p,)f%^ < rjiwrwi + \p,f%_] + c.mi + Ml- (3.44) 

From the argument of (j3.1ip and the trace theorem as well as the equation ()3.2p . 

\p,nl <\\v + iiriii < 11(1 - p)riiL + ii^iii + iiriii- 

Plugging this into (|3.44p with rj small and adding a small constant x (j3.43p . collecting terms 
and using the fact 

liprili ~ IIP/1IL> 

we obtain the uniform in e estimate 

lirilL + l/li ^ Nli + Mi. (3.45) 

We thus obtain a weak solution ^ f with the same bound (j3.45p . Moreover, we have 

e[f -f]+v- VAf - /] + Lm[f - /] = £/, [f - /]- = P^[f - /]. 

We conclude from (j3.45p . written for the difference of — /, that 

\\f-f\\l^ + \\f-f]\l < eWfWl 

< e{M\l + \\r\\l} ^0- 
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The proposition follows as e — )• 0. 

Finally from (j3.45p again we can take the limit m — )• oo of solution to 

V ■ v^r + Lmr = g, /- = p^r + r. 

By a diagonal process, there exists a weak solution / such that /™ — )• / weakly in || • 
for any fixed vtiq. It thus follows from weak semi-continuity of the norm || • H^,^^ that 



2 I l^|2 



Urn. 



+ I/I2 < 11^112 + kl2- 



Note that the limiting / so obtained satisfies jj^^-^z /(x, v)\/ ^{v)dvdx = 0. The proposition 
follows as rriQ — )• 00 and uniqueness follows from Green's identity. We remark that due to 
lack of moments control of / we cannot show /"^ — )• / strongly in L^. ■ 



4 Estimate along the Stochastic Cycles 

We define a weight function scaled with parameter 



= ^.ACW = (1 + • (4.1) 



The main purpose of this section is to prove the following: 

Proposition 4.1 Assume i3.1\} . Then the solution f to the linear Boltzmann equation h3. 2\) 
satisfies 

\\Wgf\\oO + \Wgf\oo < \\Wg g 1 1 oo + \Wg{v)r\oo. 

Moreover if g and r are continuous away from the grazing set 70, then f is continuous away 
from D. In particular, if 0, is convex then 2) = 70. 



We define 



w,{v) ^ = V27r J, (4.2) 

^e,/3,c(^)v^(^) {1 + Q^\v\^)2 

and V(x) = {v £ : n(x) • v > 0} with a probability measure da = da{x) on V(x) which is 
given by 

da = fj,{v){n{x) • v}dv. (4.3) 

We use below Definition 11.61 of stochastic cycles and iterated integral and remind the 
dependence of tk on (i, x, v,vi,V2, ■ ■ ■ , Vk-i)- We first show that the set of points in the phase 
space n^Z^Vj not reaching t = after k bounces is small when k is large. 

Lemma 4.2 For Tq > sufficiently large, there exist constants Ci,C2 > independent of 
To, such that for k = CiTq^'^, and all {t,x,v) G [0, Tq] x Q x R^, 



l{tfc(f,x,i},i}i,i}2,-,i'fc-i)>0}nj^=Jc?CTj < |-| . (4.4) 
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We also have, for /3 > 4, 



k-l 



Ei{*.+i<o<*a*.(^0(^.>n,ti^da, < |i + ^1 , (4.5) 

for all I = 1,2,..., k-l. 

Proof. Choosing < 3 sufficiently small, we further define non-grazing sets for 1 < j < k — 1 
as ^ 

= {vj £ Vj : Vj ■ n{xj) > 3} Pi {vj G Vj : \vj\ < -}. 



Clearly, by the same argument used in (j3.12p . 



daj < 



Vj-n{xj)<j 



daj + 



daj < Ci, 



(4.6) 



where C is independent of j. On the other hand, if Vj G Vj, then from the definition of 
diffusive back-time cycle (|1.35p . we have Xj — Xj+i = {tj — tj^i)vj. By (j2.4p . since \vj\ < j, 
and Vj ■ n{xj) > 3, 



(tj - tj+i) > — . 

Therefore, if tk{t,x,v,vi,V2-.-,vi:_i) > 0, then there can be at most 
S Vj for 1 < J < k — 1. We therefore have 

/ ■■■{ ^{tk>o}do'k-i} dak-2---dai 

JVi [JVk-i J 



+ 1 number of v^ 



< 5Z / ^^jZldaj 

jyi^—l 1/ {There are exactly m of fj^sVj , and k~l—m of vj^^Vj } 



< 



[a? 
E 



-1 



A;- 1 



fc— m— 1 



sup / (icTj I < sup 



do",- 



Since do" is a probability measure, Jyj ddj < 1, and 



k—m—l 



da-i 



< 



dai 



k 2 f''?^" 



But, from (^^^) < {A; - 1}™ < {/c - 1}L^^J+ , we deduce that 



^{tk>o}'^i=id(^i < 



(A; 



(4.7) 
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Now let /c - 2 = N{[^^] + 1}, so that if > 1, (gZD can be further majorized by 



r 



+1 



We choose CN,nTo^'^ ^ = ^, so that 3 = { 2CnciTo i small for Tq large and for > 3. 

Moreover, 



i- 

+ 1 ~ Cat^^Tq 



and ^'3 ° > 2, if Tq is large so that we can close our estimate. 



Finally we choose = 15. For Tq sufficiently large, '^^'3 + 1 ~ CT^^^ and k = 

^1 + 1} + 2 ~ CTq/^, and (gai) follows. 
Next, we give the proof of the first estimate in (14. Sh . The left hand side of (14. 5p is bounded 



16{ 

by 



'^{tk<o<ti}Wg{vi){vi)Yi''jZldaj 



1 f C^Ul _4 



/27r Jvi>0 (1 + ^2|^|2 



0' 



7 



where we used the change of variables : gv = u and the fact /3 > 4:. 
For the second estimate in (|4.5p similarly we have 



< / [1 + u;,(7;0(t'/>]^i^/ < J 1 + -i= / ^^l^^l—^d^; 
JVi [ V27r Ji,i>o (1 + ^2|^|2)-^ 
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where we used the fact that daj is a probabihty measure on Vj. ■ 

Denote h = Wg f and K^^^^ ■ ) = WgK{:^ •). We present an abstract iteration scheme 
which gives a unified way to study both steady and dynamic problem with diffuse boundary 
condition. Recall v = {v,v) and vi = {vi,vi) from (jl.Sp and define: 

\h'+\t,x,v)\ < lt,<oe-''^''^'\h'+\0,x-tv,v)\ 

+lt,<o [ e-''^''^^'-'^\[K^^h' + w,g]is,x-{t-s)v,v)\ds 



+lt,>o e-'^''^^'-''>\[K^X + ^o9]{s,x-{t-s)v,v)\ds (4.8) 
Jti 

+lt,>oe"'^''^^'~''^\w,r{ti,xuv)\ + —— / \Hl 

where is bounded by 

k-l 

Y.l{t,+,<o<t^}\h'^'~\0,xi - tivi,vi)\dJ:i{0) (4.9) 
1=1 

^ + Wgg\{s,xi-{ti-s)vuVi)\dT.i{s)ds (4.10) 

1=1 

k-l .ti 

+ Z] / ^{a<u}\[Kn,h^~^ + w,g\{s,xi-{ti-s)vuvi)\d^i{s)ds (4.11) 
1=1 -^'i+i 

k-l 

+ 5^1|o<tjdS[ (4.12) 



1=1 



+l{^^tu^\h'+'-''{tk,Xk,Vk-i)\dT.k-i{tk), (4.13) 
and we define 

dSz = {Ii]z]^^daj}x{w{vi)dcji}xli\-}^da, (4.14) 



Remark 4.3 For the steady case you can regard the temporal variables t, s as parameters and 
read h{t, x, v) = h{x, v). We used the notation in lll.8\) . v = {vi, • • • ,Vd; Vd+i, • • • , V3) = {v; v) 
again and v E R"' and v E R^~°' where d is the spatial dimension so that a; E $7 C R*^, for 
d= 1,2,3. 
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Lemma 4.4 There exist qq > and C > such that for all g > go, (3 > A, and for 
k = g = Cti , 

sup e'^'\\h'+\s)U 



0<s<t 
1 



< - max snp{e^^\h'-'{s)\\^}+ max \\h'+^-' {0)1 

8 0</<2fco<s<t 0<l<2k 

Wg g[s) 



+ Q 



e~^\wQ r 



C max 

l<K2fc 





h^-\s) 


/o 





{S)\oo} 

ds. 



+ sup <j e 2 

0<s<t 



(4.15) 



Furthermore for k = g = Ctt 



< e-'^^^^'-'^^Wg{vMti,xi,v)\ 



+ 



N 



+ eCp^p^N^ SUp\{vf+^r{s)\oc + CN,e,p 



sup |r(s, •, v) 



+ lsup||/i'(5)|U + e-^*sup||/i'(0)||c 

s,l 



^?sup 



Wg g{s) 



(v) 



I 



+62 sup 





h^{s) 


ds 


Jo 




2 



(4.16) 



Remark 4.5 The estimate |-^. j6p is used only in the proof of the singularity formation in 
Theorem \1.1[ 



Proof. We first prove (j4.15p and then sketch the proof of (j4.16p . 

We start with r-contribution in (j4.8l) and (14.121) . Clearly the contribution in (14. Sp is 
bounded by 

Since the exponent of is bounded by e^'^oCt-ti+i)^ from (|4.14p and Lemma[ 



1 + 



^^sWi>r{s) 



(v) 



< e 2 



_ sup <^ e 2 



Wgr(s) 



{v) 



(4.17) 



for sufficiently large > and k = g = Cti . 

We now turn to the g- contribution in (|4.8p . (j4.1Up and (|4.1ip . Rewrite 



W(,g 



{v) 



(vi)- 
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Since the exponent of dT,i{s) is bounded by e'^'^^^ the (^-contributions in (|4.8p . (|4.10p and 
()4.1ip are bounded by 



2 / e-^^''^^^-'^Wgg{s,x - {t - s)v,v)\ds 

ds 



(4.18) 



Wg{v) Jo 





W(,g{s) 


/o 


{v) 



^\ / hti+i<o<ti}U'givi){vi)U^jJdaj 
1=1 ^ 

'^{o<ti+^}Ws{vi){vi)Il]zldaj^. 



+ max 
I 



Prom 



(4.19) 



the second hne of (|4.18p is bounded by 

rt 



-!^0{tl-s) 







Wggis) 



{v) 



ds 



e 2'(*i sup <! e 2 







< g'^e 2"* sup -^6 2 

0<s<t 



u;£,5r(s 



0<s<i 

sup <5 e 2 

0<s<t 



{v) 



(v) 



With our choice k = g = Ctt and from Lemma 1121 the third and forth hne of (I4.18P are 
bounded by 



(1 + A:) 



fc-i 



/3,C 



<2{l + g){l + ^ 



/3,C (1 + 



where we have chosen sufficiently large ^0 > such that g > Qq, 



1 + 



AC 



< 2. 



(4.20) 



Therefore, the total gf-contribution is bounded by 

<;3,C (l + /)e-^* sup j 



e 2 
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Notice that the exponent m dTii{s) is bounded by e '^"^^'^ '^^ and from (|4.19p and (|4.2p and 
from Lemma 14.21 we get 

e-(-)i\h'+\0,x-tv,v)\ lt,<o 

p-v{v){t-ti) ^=-1 



+- 



+l|o<t,}|/i^+^~''(tfc,Xfc,?;fc_i)|dSfc_i(tfc) } 



fc-i 



oo / / , -^{ti+i<o<taa^« 



+/e--»* sup {e'^o^||/i^+i-'^(s)|U^ / l{o<ta^^A.-i 



0<s<t 

''^"^J ||/^^+^(0)||oo + / |l + %r max ||/.^+i-'(0)||, 

t J \<l<k 



, sup e^'''\\h^^^-^{s)\\ 

2 J o<s<t 

oo ~l~ 



jU' sup e'^%h'^^-\s)\\A , 



where we have chosen t = g = Ct^ for sufficiently large g > but fixed and used (j4.19p . 
Now we obtain an upper bound for (gS]), (gS]), (gZU]), (glTD . (gl^ . (HT^ as 



(4.21) 



Jti 



+- 



-u{v){t-h) 



E / ht,^,<o<t,}\K^,h'-\s,X,,{s),vi)\ 



+e 2 *^^(t,x,t;), 

where 74^(t, x, denotes 

= 6 2 iu;g(t;)|r(ti,xi,t;)| + -— - sup <^ e 2 



+ 



(1 + /) sup I 



e 2 ' 



g{s) 



Wg r{s) 



(v) 



{v) 



} 



1+ max ||/i^+i-'(0)||c 
J 0<l<k 



+ sup 

^ J 0<s<t 



(4.22) 



Recall from p.35p that the back-time cycle from (s, Xc\{s;t, x,v),v') denotes 

{ti,Xi,v'i), (f2, X2, U2), • ■ ■ ) i't'lh x'l' , v'l,) , ■ ■ ■ . 
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We now iterate (j4.2ip for £ — I times to get the representation for ' and then plug in 
K^gh^'^s, Xc\{s),vi) to obtain 



/R3 

-u{v'){s-si) 



(4.23) 



kwAvi,v')kw (v ,v")\h^^'^^'' {si, Xciis) - (s - si)v,v")\dsidv'dv" 



+ 



k^(vi,v')k^(v',v")\h'-^-\suX,iis) -is- si)v\v")\dsidv'dv" 



+ jj dv'dv" 



rfc-i-i)' Wg{v') 



n 1 v 



/' = 1 



<o<t;,} 



dv dv ' 



n 1 v 



dsi 1 



>0} 



+ e 2 



k^(vi,v')A^^i_i{s, X,i{s),v')dv'. (4.24) 



R3 



The total contributions of Ai^i^is in (j4.2ip are obtained via plugging ()4.23p with different 
Vs into (j4.2ip . Since J kw^{vi,v')dv' < oo, the summation of all contributions of j4£_;_i's leads 
to the bound 



</,,C 2A,.i{t) f e-^'^^'-'^e-'^'ds + Q^e-'^' max A,_,_i(t) 
Jo l<i<fc-l 



(4.25) 



fc-i 



'■{<!+i<0<t;} 



+ 



e 2 



+e--U£(f) 



1 + 



+ 



max Ai_i{t) 

0<l<k 



0<s<t 



|e^'|ti;r(s)|oo} 



0<s<t 
Q 





wg{s) 






{v) 


oo 



+e-^*|max \\h'+'-\o)\\^ + H] max sup e^^||/i^-'+i-^'(s)||oo j • 

[l<Z<2fc" [2] 0<l<ko<s<t J 

To estimate the h^~^~^ contribution, we first separate s — si < e and s — si > e. In the 
first case, we use the fact / k^gip, v)dv < +oo and by (I4.23j) to obtain the small contribution 



ee 2 max sup {e 

l<l'<k 0<<ii<s 



(^l)lloo}. 



(4.26) 



Now we treat the case of s — si > e. For any large ^ 1, we can choose a number m{N) to 
define 

v'-v"\>^. 



kmiv'i>,v") = l|„/^_^//|>J.^|j,//|<^k^p(f;,, v"). 



(4.27) 
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such that 



sup / I km (v'l, ,v")- ky, {v'l, ,v")\ dv" < ^ . 
v', Jr3 ^ 



We spUt = {kwg{v'ii,v") — km{v[i,v")} + km{v'ii-,v"), and the first difference leads to a 
small contribution in (|4.23p with k = g, 



C. 



:-fe-^^ max sup {e-f^^\\h'-^-'{s^)\\^}. 

For the remaining main contribution of km(f[,,f"), note that 

\^M',v")\<Cn- 

We may make a change of variable y = x'l, + (si — t'i,)vi' (x^/ does not depend on v'l) so that 



(4.28) 



dv' 



> e^, for s — si > e to estimate (using the notation in (|1.8|) ) 



p-i\+0\v[,\^ 

\h^-i-i (^six', + {s, - t[,)vl,v")\—-^^^^\n{4) ■ v[,\dv',dv" 



^'^ Jq J\v"\<m 



\v"\<m JV[, 

< [ e-^i+OI^^'Pdi);, / dv" 

Jn^-d J\v"\<m 



< 



' ^' {si,y,v")\dydv" 



Hence the integrand of main contribution is bounded by <£,m || — ^ — -\\2- Rearranging 



and combining KT2\\ . and KM\ we have a bound for (f03]l as 

K^X^^-\s,X,y{s),Vl) 

< e-?^[, + ^+P|Vniax sup {e^^M|/^'-'(^i)lloo} 

iV 2 J l<K2fco<si<s 



(4.29) 







t(;gr(si) 




+ Q sup 


e 2 


^^^^(•si) 




< ^ sup 










I 0<si<s 




OO- 


0<si<s 




OO- 



+ max ||/i^+i-'(0)|| 



0<K2A; 



+Ce,m,Af max 

l<Z<2fc 



/■ 


h^-\si) 


/o 


Wg{v) 



dsi 
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By plugging back (|4.25p and (|4.29p into (|4.2ip . we obtain 

I Jo Jti 

1=1 



< e-^*[£ + % + {U'] max sup {e^^||/i^-'(.)||oo} 



+e 2 '-<j ^ sup 

0<s<i 



e 2 ' 



t(7pr(s) 



(v) 



+ Q sup 

0<s<t 



e 2 



(v) 



+ max II ^^^+^-'(0) lie 
0<K2fc 



+C£^m,N max 

KK2A: 







/o 


Wgiv) 



dsi 



We then deduce our lemma by choosing t = g = Ct^/'^ and letting g large and fixed, (so are k 
and t), and then choosing e sufficiently small and N sufficiently large. It is trivial to rescale 
to the case g = 1 with a different constant depending on g. 

Now we sketch the proof of (I4.16p . Since the proof is similar to (I4.15P we just highlight 
the differences. The key is to estimate the r-contribution (j4.12p differently, using its weaker 
L\{Lf^) norm. For large > we choose m > such that 



'^\v\>m{v)wg{v)da < — . 
Prom (j4.14p and Lemma we bound (|4.12p by splitting 1 = l|^|>„i + l|^|<m: 



e 2 



k-1 



E ^o<tidak-i ■ ■ ■ dai+i 
Wg{vi)r{ti+i 



1=1 



(vi) 



6 2^'^ 



< 



k 



Nwg{v) 
^ kCjji 
Wgiv) 



^\vi\>m{vi)wg{vi)daidai_i ■ ■ ■ dai 

/ E ^0<tidcrk~i ■ ■ ■ dai+ie~^^'^+^Wg{vi)\r{ti+i)\li^^i^rnWpivi)daidai-i ■ ■ ■ dai 
1=1 

{1 + %} sup (4.30) 

J 0<s<t (v) 



\vi\<m 



\r{ti+i,xi - tb{xi,vi)vi,vi,v)\\n{xi) ■ vi\dvidai-i ■■■dai 
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We have used the daj is a probabihty measure for j > I. Using the notation (jl.Sp . ui = 
and t]j{xi,vi)vi = th{xi, -^)-^, and we have an uppder bound: 



< 



\r{ti+i,xi - ti,{xi,vi)vi,v)\\n{xi) ■ Vi\dvi 

\vi\<m 

dvi / sup\r{s,xi-ti,{xi,vi)vi,v)\\n{xi) -uiWvildvi 

\vi\<m ''\vi\<m s,v 



(4.31) 



<m / sup\r{s,xi-th{xi,vi)vi,v)\\n{xi)-ui\dui. 
Since xi — t^,{xi,vi)vi S dQ and x/ G dQ, Now apply (|2.9|) in Lemma [231 to have 



sup \r{s,xi - t]:,{xi,vi)vi,v)\\n{xi) ■ ui\dui 



snY)\r{s,y,v)\dS{y). 

an s,v 



(4.32) 



We now turn to (jOT]) . Now instead of (jl22|), combining ()i:32]) and ()i30]) yields 



e 2 '^£(t,a;,i;) 

e-?(*-*i)u;^(t;)|r(ti,2;i,t;)| 

-1 f 1 u'er(s) 

.4 Hiv^^^p 



w}g(u) 



(v) 



+ Cn 



sup |r(s, •, v)\ 



+C/3,c(l + /)sup 



5(s) 



+ e--* sup ||/i'(0)||oo + C^,c sup ||/i'(s 



''^^'-'''>w,{v)\r{h,xuv)\+D. 



(4.33) 



We now plug (j4.33p back into ()4.24p . the estimate for (j4.23p . for each £. Now we obtain the 
upper bound for the last term of (|4.24p using the estimate of ^£_i_/(s, Xci{s),v') from (|4.33p 
as 



e 2 ' 



k'e'f''^w,{v')\r{t[,x[,v')\dv' + e-'^'CKD 



R3 



where 



is the exit time of the point {Xc\{s),v') and 

x'l = Xc\{s) - th{Xc\{s),v)v, k' = k^, (Vc\{s),v'). 



(4.34) 
(4.35) 
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We further plug this bound back into (j4.23p and then back into (j4.2ip to isolate only the 
contribution of Following exactly the same steps, it suffices to control 

U,<o /*e--('^)(*-«)e-^^ / k'e'^''^w,{v')\r{t[,x[,v')\dv'ds (4.36) 



+lt,>o e-'^^'')^*-^) j k'^^e-''^w,{v')\r{t[,x[,v')\dv'ds 



+ , , . X 

Woiv) 



+ / l{o<t,+i} / K,,(i^^'^Wi,{v)\r{t[,x\,v')\ \di:i{s)ds 



+e-^(*-*^)u;,(t;)|r(ti, XI, + I? + { [1 + + ^} ^' 

by Lemma 1121 Note |[1 + ^^]^ + ^^j < 1. We shall use a sequence of approximations 
to estimate the above integrals. Let k!^^l^^,_y |„/_yj(5)|<m — ^wg- We again split for 

m large 

/ Ik' _ V" I < — 

j \^we ^"-el - jV' 

We therefore bound 

+lt,>o I' e-'^^^^'-'^ I K/-^''^w,{v'Mt[,x[,v')\ds 



^^■^■'•^■sup|t^,r(5)|oo. (4.37) 



Since k™^ is bounded, and Wg{v') <g '^^ further split > N and < N. Since 

i|t;|>Af (^)~^ ^ ^ constant of CN,s,k the main integrals in (|4.37p are bounded by: 

lti<o / / {v')'^\r{t'i,x'i,v')\dv'ds + lt^>o / (u')^|r(t'i, x'^, t>')|df'(is 

Jo JUm Jti JUm 

+ E| /'i{*..i<o<*.} / {vr\r{A,x',y)\ 

+ r l{o<t,+i} / (^;nr(t;,x;,z;')l|dSKs)ds 
+^sup(t;)^+V(s,a;,t;)|, (4.38) 

-1' s,v,x 
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where 



v'eUm iff \v' -Vci(s)\> —, and \v' - Vci(s)\ < m, and \v'\ < N. (4.39) 

m 

Lastly, for any e > 0, we further spht the time intervals of the main integrals in (j4.38p to 
obtain: 



rt-e 



lti<0 



+ 



{v')'^\r{t'i,x\,v')\dv'ds + U^^o 



rt-e 



{v')^\r{t'i, Xi,v')\dv' ds 



Iti+e JUn 



k-l 



ti-e 



+^Cg^k,i3 sup(v)'^|r(s. 



+ 



X.V]\. 



1+ 
ti-e 



l{o<t,+,} / {vYHt'^,x'^,v')\\dT.i{s)dv'ds 



(4.40) 



s,x,v 



We now are ready to use change of variables to estimate the main f'-integrals in (|4.40p . Recall 
that, by the definition (|4.35p of x\ and Definition ()1.6p . X(.\{s) reaches if and only if s at 
these t,ti,t2, ■■■ti^i. From our splitting of time intervals, there exists > such that 

dist(Xci(s),(9il) > Ce > 0, 

in the integrals in (j4.40p . uniformly in Hat. We now repeat the change of variables x'l — ?■ ||7y 
as (OB and (02D . but using (|2J0]1 instead of (gj]). We finally bound (ICTD by 



C 



N,e,k,e 



sup |r(s, •, v) 



+ £Cg^k,p,N s\xp{v)^\r{s,x,v)\. 



(4.41) 



Collecting and combining (|4.33p . (|4.36p . (|4.37p . (|4.40p and (|4.4ip . we conclude the r contri- 
bution in (j4.2ip is bounded by 



+CN,e,k,s 



snp\r{s,-,v)\ 



We deduce (|4.16p by recalhng D defined in (|4.33p . ■ 

After preparing the above tools we return to the stationary problem to give the proof of 
Proposition 14.11 

Proof of Proposition 14.11 We use the exactly same approximation ()3.5p to establish the 
proposition and follow the same steps in the proof of Proposition 13.11 We denote /i^"*"^ = 
Wg /^^^) and rewrite (I3.5p as 



Wgg, 



(4.42) 



Jn(x)-v'>0 



h^{t, x, v')wg{v')da + Wg r. 
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Step 1: We take ^ — )• oo in L°°. Upon integrating over the characteristic hnes ^ = v, and 



dv 
dt 



and using the boundary condition repeatedly, we obtain (by replacing 1 with 1 — 7 



and V with v + e) that the abstract iteration (j4.8p is valid for stationary h^^'^{s,x,v 
h^^^{x,v), g{s,x,v) = g{x,v) and r{s,x,v) = r{x,v). Therefore, for £ > 2k, by Lemma 
we get (choosing k = g = Ct^/^ large) 



^^"^^1100 < - max {II /i^ '||oo} + e '2°* max \\h'^^^~ 



l<l<2k 

+Q 





Wg g 




\wq r\oo + 


{v) 






00- 



0<K2A: 



+ Ck max 

l<l<2k 



Then, absorbing e 2°*||/i^+i||^ in the left hand side, we have 



< - max {\\h' 

4 l<«<2fc 

+Ck max 

l<l<2k 



'+1-1 \ 



{wq r 00 + 


g 




{v) 






CO- 



I- 1 



Now this is valid for all i > 2k. By induction on £, we can iterate such bound for i+2, ....i+2k 
to obtain 



||/i'+i|oo < 7 max {\\h'+'-'\\oo} + Ck 

4 l<K2fc 



< I max {||/i'-^+*-'||oo} + 2Cfc 

4 l<K2fc 





wg g 




{v) 





\Wg r\oo + 


Wg g 




{v) 






00 



+ max 

-2k<l<2k 



+ max 

-2k<l<2k 



< \ max {||/i^-'||oo} + (2i + l)Cfc 

4 1<«<2A: 

We now take a maximum over 1 < i < 2/c to get 



\Wg r\oo + 


wg g 


{v) 





+ max 

-2k<l<2k 









2 



max ||/i^^"^ 'lloo ^ ~ max \\h 

l<l<2k 4 l<l<2k 

+ {Ak + l)Ck 



Wg r|oo + 



Wg g 



(v) 



+ max 

-2k<l<2k 



This implies that, from Lemma 14.41 for all i > 2k, 



(4.43) 



max \\n ■ \\oo 
l<K2fc 



max \\h lloo + \wg r\oo + sup 

l<'<2fc 0<s<t 



Wg g 



{v) 



+ max / 2 

1<K£ 



(4.44) 



Now in order to control maxi<K2fc ||^'| 
to obtain, by Lemma 14.21 



we can use (j4.8[) repeatedly for h'^^ — ^ h"^^ ^ 



max ||/i'||oo <fe I'W^, r|oo + 

l<K2fc rs.. I t/ I 



(4.45) 
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We therefore conclude that, from (j4.44p and (j4.45p . 

,0| 



max ll/i^"*""^" 

l<l<2k 



oo ^fe 



0<s<t 



Wq g 



{v) 



+ max llf II2 

KKl 



\h Woo + \Wq r\oo + sup 
But maxKKoo ll/'lb is bounded by step 1 in the proof of Lemma l3.31 Furthermore for /? > 4, 



is bounded by \\wg ■ ||oo and \r\2 is bounded by \wQ{v)r\cyo. Hence we have 



max ^lloo <h 

l<K2fc 



+ \wq r\oo + sup 



0<s<t 







{v) 


00- 



Therefore there exists a hmit (unique) solution ^ h = Wg f € L°°. Furthermore, h satisfies 
gS]) with /I'^+i = h. 

By subtracting h^~^^ — h in (j4.8p with r = (7 = 0, we obtain from Lemma 14.41 for h^~^^ — h, 
(choosing k = g = Ct^/^ large) 



\h'+^-h\\ 



< 



< 



■ max 1 1/1 

0<K2A; 



T II -!^k II 7 / 

n 00 + e 2 max \\h 



0</<2fc 



+C7fc max ||f-'-/||2 

l</<2fe 



1 



max 



4 l<l<2k 

where e"^''||/i^+^ - 00 is absorbed in the left hand side 



+ C7, max ||f-'-/||2, 



< 



< 



1 



max \\\h! 



4^/2fc l<l<2k 
Ck{\wgr\oo+ supo<^<( 



max - /II2 

-2fc<K2fc"'' 



Wgg 



4£/2fc 



max - /II2 

-2fc<z<2fc"-' -"'^ 



From step 1 of Lemma [3T3l we deduce ' — /i||oo — for ^ large. 

Step 2: Now we let j — )• 00. We take to be the solution to ()3.7p and integrate along 
dx _ „, dv _ Q repeatedly. We establish (j4.8p for /i^ = /i-^' (we may replace (1 — i) by 1 and 



dt 



dt 



£ = to preserve inequality) and / = 0. Lemma 14.41 implies 



\wgr\oo + 



Wgg 



so that 



\W\\oo<Ck 



Wgr\oo + 



w 



o9 



(v) 



(v) 

+ Ck\\fh. 



+ Cik)\\p 



Since ||/-'||2 is bounded, this implies that ||/i-'||oo is uniformly bounded and we obtain a 
(unique) solution h = wf G L°^. Taking the difference, we have 

£[h^ -h]+v V,^[h^ -h] + v[y -h]= K^[y - /i], 

y_ — h- = , , ^ f [h^ — h]{t,x,v')wg{v')da{v') 



y (t, X, v')wg{v')da{v'). 



Wg{v) Jn{x)-v'>0 

1 1 
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We regard — j ^ ^^^^ f^f^^y,^,^Qh{t,x,v')wg{v'){n{x) ■ v')dv' = r. So Lemma implies that 



\W - h\\oo < - h\\oo} + ^W\oo + Ck [\\P - fh] , 

which goes to zero as j to oo. 

We obtained a solution /i*^ = Wgf^ to (j3.3p . By integrating over the trajectory, (14.8 
is valid for replaced by /i^ so that from Lemma ST 



and hence 



loo < l{\\h' Woo} + e-^'Wh'W^ + k 



\h'\\oo<k 



\wr\oo + 


wg 




(v) 






OO- 



+cik)\\n\2, 



\wr\oo + 


wg 




(v) 






OO- 



which implies that, from Proposition 13.11 that ||/i'^||oo is uniformly bounded and we obtain 
h = wf solution to the linear equation. 
Now we have 



eh'^ + v ■ Vxlh'^ -h] + vLm[h!^ - h] 

ht - /i_ 



0, 



— h]{t,x,v')wg{v')da, 



so that from Lemma HI 



'^giv) Jnix)-v'>0 

\\h' -h\\oo<ke\\h'\\oo + \\r -fh, 

which goes to zero. For hard potential kernel, we recall that in previous section we have 
constructed an approximating sequence Z™" to the equation 



9, 



with uniform bound in and a limit — )■ / weakly in || • see Step 2 in the proof of 
Proposition EIT] Moreover, we obtain from ()4.4|) that ||'Wg/'"||oo is uniformly bounded and so 
is Wgf . Note that 

V ■ v.[r - /] + L^ir - /] = [-L + L„]/, [r - /]- = p,[r - /]• 

It follows from Proposition 13.11 and the boundedness of Wgf that 

iir-/ii.<ii[-i^+L^^]/ii2^o. 

Now, to show 11/™ - /lloo ^ 0, we apply with g = (-L + 

{Lm - L}f 



l|{/'"-/}||oo<||r-/||2 



{v) 



0. 



Since Wgf £ L°° , the second term goes to zero. 

In the iteration scheme is continuous away from D and hence / is continuous away 
from ■ 

Remark 4.6 Our construction fails to imply the Fg = ^ + yjjlfg > 0. This can only he shown 
by the dynamical asymptotic stability discussed in Section 7. 
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5 Well-posedness, Continuity and Fourier Law 



Proof of the Theorem 11.11 

Wellposedness. We consider the following iterative sequence 



(5.1) 



1 , N-fj- 



J n{x)-v>0 



f^^/Ji{n{x) ■ v)dv + 



with /o = 0. 

Note jT{f\ f)y/Jl = and from /„(^).„<o ' ^')^^ = L{xyv<o t^{n{x) ■ v)dv = 1, 



7_ I Jn(x)-v>0 



f^{n{x) ■ v)dv + ^^^^ J> d7 = 0. 



Since \wg ^^^ |oo ^ (5, we apply Proposition 14. 1 1 to get 



'>J ||00 I "JQj |00 



{v) 



Since 



(v) 



^ Ike/ II we deduce 



so that for 6 small, 



)./ lloO I "^pj loo 



Upon taking differences, we have 



= nf-f-\f)+nf-\f-f-') 



[f-f-\n{x)-v)dv + ^, 



J n(j;)-f >0 

and by Proposition 14.11 again for f^~^^ — 

\\w,[f^^ - f]\\oo + \w,[f^^ - /]ioo < 5{\\w,[f - f-^]u + kjr - r-']ioo}. 

Hence is Cauchy in L°° and we construct our solution by taking the limit — )• fg. 
Uniqueness follows in the standard way. Moreover due to Theorem 2 and Theorem 3 in |26] 
is continuous away from T) and so is fs- Moreover, if is convex, then J) = 70. 

Formation of singularities. Now we prove the formation of singularity. Recall that li^loo ^ 1 
and the wall Maxwellian is defined as 



1 



27r[l + 



■ exp 



2[1 + 5^{x)] 



43 



1 I I 

while the global Maxwellian is = 7^ e 2~. Then 

fis{x,v) - fi{v) - 2[— l]i^{v)6^{x) 

s:2i Q|2 



< 5^|^|^[l + |Hlexp 



2{l + mc 



(5.2) 



For any non-convex domain $7 C R"' for d = 2,3, there exists at least one {xq, v) G 7^ with 
/ 0. Suppose h = w^f satisfies (|5.ip with f^^^ = = f- We prove that h is discontinuous 
at {xo,v) by contradiction argument. Assume h is continuous at {xo,v) where the magnitude 
of V will be chosen later. By definition 



, V . V 
Xi = Xo- ti,{xo, 7-r)737. 

\v\ \v\ 

We choose to be a continuous function such that 



(5.3) 



and mainly zero elsewhere. Then by the continuity assumption we know that the quantities 
I and II defined below must coincide: I = II. I is obtained by evaluating h{xQ,v) through 
the boundary condition 



h{xo,v) 



1 



Ws{v) Jn{xoyv'>0 
+ Ws{v) ^ \{xo,v)- 



h{xo, v')ws{v')d(j 



7l(xo)-f '>0 



h{xQ, v')w{v')da 



(5.4) 
(5.5) 
(5.6) 



On the other hand, the existence of xi allow us to evaluate h{xQ,v) = II along the trajectory 
and has the expression (14. Sp with H as in (14. 9p - (I4.13P with 

r = — — , h =h, g = r{ — ,—), 

yTX W W 



which are independent of time. Since from (j5.2 



sup |r(-, v)\ 

V 

Wg g 



r\2 < 0\V\2, 



1 



{v) 



< 



{v) 



<\\^ph\\l<S'ml<6' 



2 < S'\^\t,+5m, 
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where we have employed Propositions 13. l l and 14. ll for the last two estimates. Due to (j5.2p we 
can use the estimate (|4.16p for II : 



\h{xo,v)\ 



Wg(v)\r[xi,v) 



\oo + Q 



Wg g 



(v) 



sup \r{-,v^ 

V 

t 



+ C 



, ds 







(5.7) 



= e-''^^^^'-'''>w,{v)\r{xi,v)\+B. 

Now we will show that, for a suitable choice of I — II ^ 0, which is a contradiction so 
that h = Wgf has a discontinuity at (xq,v). 
Rewrite 

I-II > 

v 

-{\^\ + \&\+B}. 
For large by ([52]), t > ti, e"''('')(*~*i) < 1, we have 



(5. 



.-u(v)(t-u).^^^.. fs{xi,v)-ti{v) 



C5' 



I |2 I |2 

> w,{v)[^-^-l]^6\^\oo-C6\ 
Recall that ||/i||oo ^ 'JI'J^loo ^ <5 to conclude that 

I rai .1(1531)1 

I 1 2 

By (|5.9p . we can find |fo| > so that c^q > sufficiently large so that 



(5.9) 



\vo\ 



1] \//U(fo) - '*U'g(t>o) VlAvoj > Cva > 0. 



On the other hand in ()5.7p and ()5.9p . we choose 6 sufficiently small, and p sufficiently large, 
then N sufficiently large to get, then e sufficiently small, such that for v = vq, 

I_II >^S-C6\^\2. 



Since -d is almost zero except for xq and xi, we can can make |'i?|2 arbitrarily small. In 



particular, there is a continuous function -i? such that ^6 — C6\'d\2 > > 0. Hence 



I — II > and this is a contradiction. 
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Proof of Theorem 11.21 

We now prove the (5— expansion. Since ^s{x) is analytic with respect to 5 we have 

Hs{x,v) = 11 + 5^1+ 5^ ^2 -\ • (5.10) 

Note that ^ Pi{v)e 2 , where Pi{v) is some polynomial. Further we seek a formal 

expansion 

Plugging this into the equation, 

V ■ V,[5/i + d^f2 + ...]+ L[5h + 5^2 + . . . ] 

= r(5/i + ,52 /2 + ...,<5/i + 52/2 + ...), 

[<5/i + 6^ f2 + ■■■]- = Py[Sfi + <5'/2 + . . . ] 

+ \^!^l±f^l±llA I + sf, + 6'h + ... ] . v)dv. 

Jn{x)-v>0 

We compare the coefficients of power of 6 to get an equation for /j for i = l,...,m — 1, 
(assuming /q = 0) 

i-l 

vV.h + Lfi = Y,r(fj,n.j), (5.11) 
i=i 

i— 1 „ 

/ik = Pyh + Y.-^ h^,^{n{x)-v)dv + 

j^l V/^ Jn{x)-v>0 



Note that, from the (5— expansion, 

/ {^^ {x,v) — ^{v))\n{x) ■ v\dv = 5^"^* / /ii|ra(x) • v|(it; = 0. 

Since 

i— 1 i— 1 

iR3 ^.^-^ Jn{x)-v<0 ^ ^/f^ Jn{x)-v>0 y/^ 

by applying Proposition 14.11 repeatedly, we can construct /i, /2, /m-i inductively so that 
for ^ C < i, 

\\Wgfi\\oO + \'Wgfi\oo <ml- (5.12) 

In particular fi satisfies (jl.lOp . Now we define the remainder satisfying 

fs = Sh + --- + 6"'-'fm-i + 6^fi, 
and we obtain the equation for the remainder 

+ r(/i, /i + 5/2 + . . . + ,5— + <5'"r(/4, /;^) + /, 

fLU- = P,fm + ^^^l fi^{n{x)-v)dv + r'. 

Jn(x)-v>0 
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where 



g 



2m— 1 m—1 iiL „ 

i=m 7 = 1 i=l V Jn 



i=m j=l 

Since by ([51^ 



' n(x)-7;>0 



fm-i^/Jiinix) ■ v)dv. 



Wg g" 



{v) 



+ \Wg r 



loo 1) 



we can apply Proposition 14. II to deduce that 



1"^^ /ml loo "I" \ Wg fm\oo 



<m S\\Wg /4||oo + S"^\\Wg f^Wl^ + 6\Wg | oo + \\w g \\oo + \w g / | oo • 

We therefore conclude that \\wg fmWoo + \ wg fm\oo ^ 1 and the expansion is valid. ■ 

Proof of Theorem 11.51 We now consider a slab — ^ < 2; < We consider the stationary 
solution fs G satisfying 



vAfs + Lfs 
fs{x,v) 



TifsJs), 



(5.13) 



J n(a;)-f >0 



fsy/Jl{n{x) ■ v)dv, 



for X = —\ 01 X = From Theorem 11.11 \\wg fs\\oo ^ <J- We claim that dxfs{x,v) G 
i + e, +i — e] X R^) for any small e > 0. In fact, we multiply (I5.13|) by a 1-dimensional 
spatial smooth cutoff function x{x)i (x = ii^ar x = — ^ and x = +i) so that 

vidMs] + Lixfs] = Tixfs, fs) - vifsx'- 
We take one spatial derivative (jS.lSp to get 

vAxixfs] + Ldx[xfs] = T{dx[xfs],fs) + r(x/s, dxfs) - vix'dxfs - vix" fs 
= r{dx[xfs]Js) + Tixfs, dxfs) + x'iLfs - T{fs, fs)) - vix"fs 

= r{dx[xfs]Js) + r{fs,dx[xfs]) 

-r(/„ x'fs) + x'{Lfs - r(/„ /,)) - vix"fs, 

where we have used (j5.13p to replace —vidxfs = Lfg — T{fs, fg). Letting Z = xdxfs, then we 
have 

viZx + LZ = r{zjs) + r{fs,z) 

-nfs, x'fs) + x'iLfs - r(/„ /,)) - v,x"fs. 

Note that -r(/„ x'fs) + x'{Lfs - r(/„ /,)) - vix"fs G Multiply by Z and use Green's 
identity to have (no boundary contribution) 

\\{I-P)Z\\l < d{ \\Z\\l + 1 }. 

Then we repeat the proof of Lemma 13.41 but replacing (/)a in (I3.37P with the solution of 
— A(/)a = a with (pa = on dQ. Since the boundary condition of Z is Z^ = 0, Lemma 13.41 is 
valid so that 

WxdxfsWl < \\Z\\l < +00, 
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and docfs is locally in L^^^i^, 1). 

Recall that ^ 

'^s{x) = -n — r-j- vFs{x,v)dv. 

Now we are going to show 

Us{x) = 0, for any X e [-^,+^]. 

For any smooth test function ^ G C°°([— 5, +|]), due to the weak formulation for Fg, we have 
for the first component of pgU that 

/ ^ V {x){psUs)idx 

[H , r 

= ip (x) viFs{x,v)dvdx 
J -I iR3 

= V'(^) j^^ v,Fs{\) - ^{~) J^^ v,F,i~) - jy i^ix) J^^ vAFsix, v) 
= -['2p{x)f Q{Fs,Fs){x,v)dvdx 



2 

0, 



where we used the boundary condition as well as the orthogonality of Q to the collision 
invariants. Hence dx[{psUs)i{x)] = in the distribution sense and hence pgUg G 1^^'°°. 
Therefore pgUg is continuous up to the boundary and lim^^_ 1 psUs{x) = psUs{—^) and 

lim^^^i PsUs{x) = psUs{+^). Moreover for all i/j G C°°{[-1, +^]) 

^ lk'{x){psUs)l = {psUs)l{^)^{^) - {psUs)l{-^)lk{-^), 

to conclude 

1 1 

{psUs)i{x) = 0, for all x G ["2' +2^' 

On the other hand 

Ps{x) = / fi{v) + ^/JIfs{x,v)dv > / p{v)dv - / ^{v)w~'^{v)dv X \\wfs\\oo > 0, 

for small ||tf/s||oo so that ps is positive. Thus {us)i{x) = 0. The components {us)i for i = 2,2, 
vanish by symmetry: 

{us)i{x) = I ViFs{x,v)dv 

Ps[x) Jr3 

= — / ViFs{x,v)dv ^ / ViFs{x,v)dv 

Ps{x) Jvi>0 Psix) J^.^Q 



Ps{x) Jvi>0 ' Psix) Jv^>o 



ViFs{x,v)dv --^ / ViFs{x,v)dv, (5-14) 
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where {v)i = —Vi and vj = vj for j ^ i, and we have defined Fs{x,v) = Fs{x,v). Then F 
solves 

vA,Fs = Q{Fs,F,) , Fsix,v) = fi^{x,v) [ Fs{x,v){n{x) ■v)dv, 

where x E 50 and n{x) • u < 0. By the uniqueness, we conclude Fs{x^ v) = Fs{x, v) and hence 
Fs(x,v) = Fs{x,v) almost everywhere. Therefore it follows that (j5.14p vanishes and {us)i = 
almost everywhere for i = 2,3. 

Recall the definition of temperature Og in (|1.23p . Since dxfs G | + e, + ^ — e] x R^), 

we have dx9s £ + +f — e] x R^) so that in the sense of distribution, 

with K' = K > 0. Since Ug = 0, by Green's identity, 

dxll; v\v\^fs \ = dxQs = 0. 



/R3 

By the Fourier law (jl.ip we have then, in the sense of distribution 

dxAK{9s)} = 0. 
Therefore, in the sense of distributions 

K^Os) = Ax + B. 

But if 6*, = / blVs^/i" G then K(6l^) G and this implies that i^(6's) is in W^^.oo 
so that K{9s) is continuous up to the boundary. But if = k > we thus deduce that 9^ 

2' ^2' 



itself is continuous up to the boundary x = —•k^+h. We then can rewrite 



By the main theorem, fg = dfi + 0(5^), so that 

9s{x) = 9o + 69i{x) + 0{5^), 



and we know that ^0 = 5/ 1^1^/^^ is constant in our construction. Then we have for all 
K{9o + S9i + 0{d^)) = K{9o) + 5K'{9q)9i + ©(.J^). (5.15) 



'2 — — ~'~2' 



We therefore have 

A = 5K'{9,)[9,{+\)-9r{-\)] + 0{6^), 

B = Ki9o) + 6K'{9o) ^'^^^^ ^ ^'^"^^ + 0(6^). 
Hence the first order expansion 9i{x) must to be a linear function: 

Here 9i = / bP/i and from (fOTT) /i satisfies ([LM]) . ■ 
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6 Decay 

The main purpose of this section is to prove the following: 
Proposition 6.1 Suppose that for all t > 

g{t,x,v)y/JIdv = r^d'y = 0. (6.1) 



Then there exists a unique solution to the problem 

dtf + v-V,f + Lf = g, /(0) = /o, mOxR3xR+, (6.2) 

with 

f-=P^f + r, on7_xR+, (6.3) 

such that for all t 0, 

// f{t,x,v),/jldxdv = 0. 
J Jnxn^ 

Moreover, there is X > such that 



(6.4) 



Wfml + \fml < e-''[\\fo\\l + 1^ e'^g{s)\\lds + e'y{s)\lds}. 

Lemma 6.2 Assume that g and r satisfy ^6.1\) and f satisfies ^6. S\) , lid. 3\) and (6^. Then 
there exists a function G{t) such that, for all < s < t, 

t\\Pf{r)\\l < G{t)-G{s)+ f\\g{r)\\l + \riT)\l 

J s J s 

+ r||(I-P)/(r)||2+ f\{l-P,)f{T)\l^. 

J s J s 



Moreover, G{t) < \\f{t)\\l 

Proof. The key of the proof is to use the same choices of test functions (allowing extra 
dependence on time) of ()3.16p . ()3.23p . ()3.30p and (I3.37P and estimate the new contribution 
Is I /cxR3 ^t^/ time dependent weak formulation 

* [ i^fdj - f I ijfd^ - f 11 ^ • v.^/ 

s J-fj^ Js J'-f^ Js J JQxII^ 

= - [[ i'fit) + If i^fis) + f II -m^ - p)/ + H 

+ f II dti^f- (6.5) 
Js J JnxR^ 

We note that, with such choices 

G{t) = -ll 4^f{t), \G{t)\<\\fml 

J JnxR3 
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Without loss of generality we give the proof for s = 0. 

Remark: We note that (j6.ip . (j6.2p . (j6.3p and (|6.4p are all invariant under a standard t- 
mollification for all i > 0. The estimates in step 1 to step 3 are obtained via a t-mollification 
so that all the functions are smooth in t. For the notational simplicity we do not write ex- 
plicitly the parameter of the regularization. 

Step 1. Estimate of V x^^dta = Vxdt4>a in (|6.5p 

In the weak formulation (with time integration over [i, t + e]), if we choose the test function 
ip = 4>y/JI with dependent only of x, then we get (note that Lf and g against <l){x)^/J2 
are zero) 

/• rt+e p rt+e r 

V2^ / [a{t + e)- a(t)](/.(x) = 27r\/2^ / / (6 • V^.)0(a;) + / / ^0^, 
Jn Jt Jn Jt Jj^ 

where /j^a fi{v)dv = J-^3{vi)'^fi{v)dv = 2'K\phi and we have used the splitting ()3.19p 

and (|3.20p . Taking difference quotient, we obtain, for almost all t. 



/ <^dta = V2^ f {b- Vx)(t> + ^ / rcp^. 

u S~2 t/ t/ ^ 



Notice that, for (p = 1, from (|6.ip the right hand side of the above equation is zero so that, 
for all t>0, 

dta{t)dx = 0. 



in 

On the other hand, for all <j){x) € = H^{Q), we have, by the trace theorem \(j)\2 < ||(/>||j|/i, 



n 



(j){x)dtadx 



< 



< 



kl2|0|2 + ||6|k 



2\m\m 



{ m)\\2 + \r\2] 



H 



1 . 



Therefore we conclude that, for all t > 0, 

\\dta{t)\\^H^y < ||6(t)||2 + |r|2, 
where (H^)* = (i/^(r2))* is the dual space of H^{Q) with respect to the dual pair 

{A,B) = [ A{x)B{x)dx, 
Jn 

for ^ e and 5 G {H^f . 

By the standard elliptic theory, we can solve the Poisson equation with the Neumann 
boundary condition 

- = dta{t) , 



= 0, 

an 



dn 

1; 



with the crucial condition dta{t, x)dx = for all t > 0. Notice that = ~ dt4>a 



where 4)a is defined in (13.37p . Moreover we have 

||Va;(9t(/'a||2 = \\^'^^dta{t)\\Hi = W^aWm 



< 



\\dtam^jjry<\mh + \r\2. 
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Therefore we conclude, for almost all t > 0, 

\\V^Mam2<\m\\2 + \r\2. (6.6) 



Step 2. Estimate of Va:/^~^dtV = V^dtct^l in ([63]) 

\v\'^ 

In (|6.5p we choose a test function ip = (f){x)vi-y/JI. Since fi{v) = -^e"^", J ViVjfi{v)dv 
/ ^)fj,{v)dv = 2-7TV^6ij and we get 

27r\/2^ / [6i(t + e) -6i(t)](/. 
Jn 

= - / f(t)Vi^ + 27r^/2^ / / di(t)[a + c] 

Jt J-y Jt Jn 

+ / // y^v,Vi^d,(t>{l-V)f + II (t>v,g^. 

Taking difference quotient, we obtain 



n 

-1 



^ ,^ , f{t)vi^^/JI+ / ma{t)+c{t)] 
Ztt\/2tt j^ jq 



1 



y^vjVi^dj(l){I-P)f{t)+ [[ ^Vig{t)^}. 
3 JJnxRS J 



27r-\/27r I- J JqxR^ ^ 

For fixed t > 0, we choose (p = solving 

-A$^ = a<6,(t), ^ilan = 0. 

Notice that <I>^ = — A^^5t6i = where (pl is defined in (|3.24p . The boundary terms vanish 
because of the Dirichlet boundary condition on Then we have, for t >0, 

[ \V,A-^dMt)\^dx = [ iV^^ifdx = - [ A^^dx 
Jn Jn Jn 

< ^{l|v.$illi + ll'J>^lli} + lla(t)ll^ + l|c(t)||i 

+Il(i-p)/(t)i + ll5(t)iii 

< e\\v,m\i + \\ami + \\cmi 

+ \\il-P)fit)\\l + \\git)\\l 
where we have used the Poincare inequality. Hence, for all t > 

Wv.Mmh = \\v,A-'dMt)\\2 

< \\a{t)h + ||c(t)||2 + 11(1 - P)/(i)l|2 + imh- (6.7) 



Step 3. Estimate of V^A'^dtc = V^dtcpc in dSS 
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In the weak formulation, we choose a test function (j){x){^-^ Since J-^j 

-)dv = 0, /j^3 iJ,{v)viVj{ ^^^ = 27r\/27r(5jj and J^-^ fi{v){^-^ ^)'^dv = Zixyf^ we get 



2 



Jvl Jn 

/t+e r- i-t+e r- 1 |2 o 

+ [[ (I - P)/(M! _ . V.)0 

and taking difference quotient, we obtain 

4>[x)dtc(t^ x)dx 

- [ b{t) ■ V^cl) [{—- -)Mf{t) 



1 ff ,^ „..JwP 3. 



37rV27r J jQxRS ^ 2 



Iff . , . JvP 3. 



37rV27r J JnxR3 ^ ^ 

For fixed t > 0, we define a test function = <I>c where </)c is defined in (I3.16p . The boundary 
terms vanish because of the Dirichlet boundary condition on <I>c. Then we have, for t > 0, 

-A$, = dtc{t), = 0. 

Notice that $c = -A~^dtc{t) = dt<t)c{t) in ([3T6]) . We fohow the same procedure of estimates 
Va;A~^(9ta and Va;A~^(9t5 to have 



Jn 



^c{x)dtc{t, x)dx 

n 

2 I ll/F. I|2t I ll;,/'+M|2 , ii/t T^^ -P/'4-M|2 



< e{||V.$,||^ + ||<I>,||^} + + 11(1 - P)/(t)||^ + 

< e\\V,^,\\l + \\b{t)\\l + \\{I-P)f{t)\\l + \\gml 

where we have used the Poincare inequahty. Finally we have, for all t > 0, 



< ||6(t)||2 + ||(I-P)/(t)||2 + ||5WI|2. (6. 



Step 4- Estimate of a, b, c contributions in (16.5P 
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To estimate c contribution in (j6.5p . we plug (j3.16p into (j6.5p to have from (j3.20p 

f ff {\V\^ - ^3c)v^^^t^^ct)J 

Jo J JnxRS 

= XI / / / (1^1^ ~ f3c)viVjti{v)dtdi(j)cbj 
j^-^ Jo J JnxR3 

Jo J Jqxrs 

The second Hne has non-zero contribution only for j = i which leads to zero by the definition 
of in (|3.2ip . We thus have from (|6.8p . for e small, 

r //' {\v\^ - (3,)v,^dtdi(i)j 

JJuxH^ 



< 



< 



^ {||6||2 + ||(I-P)/||2 + ||5||2}||(I-P)/||2 

11^112 + ^*{ll(I-P)/lli + ll9lli}. (6.9) 



Combining with (j3.22p in Lemma 13.41 we conclude, for e small, 

c{s)\\lds < G{t)-G{0) (6.10) 

+ ^* {ll(I - P)/(^)ll^ + Us)g + 1(1 - P,)f{s)\l+ + \r{s)\l + e\m\\l}ds. 
To estimate b contribution in (|6.5p . we plug (|3.23p into (j6.5p to have from (j3.20p 




JJQxR^ 

= f II (^.'-/3.)MWjM!_|}e (6.11) 
Jo JJnx-R? ^ ^ 

+ 1 II {vf-Pb)VJ^dtd,4{i-p)f, 

Jo JJnxR^ 

where we used (I3.28P to remove the a contribution. We thus have from (16. 7p . 

f II {vf - /3b)VJ^dtd,cl>if 
'0 JJnxR^ 



< 



< 



^ {llalb + llclb + 11(1 - P)/||2 + lbl|2}{||c||2 + 11(1 - P)/||2} 

l^{\\{I-P)fg + \\c\\l + \\gg+e\\a\\l}. (6.12) 
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Next we plug (|3.30p into (j6.5p and we have from (j6? 

n\v\^ViVj^dtdj(j)lf 
JxR3 

= f [ \vWvj^dtdjcpi{i-p)f 

Jo JQxR3 



< 



[ {||a||2 + ||c||2 + ||(I-P)/||2 + N|2}||(I-P)/||2 

Jo 

< ^* { 11(1 - P)f\\l + Ml + e[\\ag + \\c\\l] }. (6.13) 
Combining this with (|3.36p in Lemma 13.41 we conclude from (j6.5p that 

f \\b{s)gds < G{t) - G{0) (6.14) 
Jo 

+ { 11(1 - P)f{s)\\l + ||5(.)||2 + 1(1 - P,)f{s)\l^ + |r(.)|2 + e\\a\\l + ||ci(l + e)}ds. 

Finally in order to estimate a contribution in (|6.5p we plug p.37p for into ()6.5p . We 
estimate 

(b|'-/3a)^^^M9^<^a/ (6.15) 

Jo Jqxr^ 

= f f {\V\^ -(3a){VifM^ahi 

Jo Jqxrs 



+ 



f [ {\v\^-f3a)VifldtdiMi-'P)f 

Jo Jnxn.^ 



< 







{\m\\2,n + \r\2}{\\bh,n + \\il-P)fh} 
Combining this with Lemma 13.41 we conclude 

f\Hs)\\lds < G(t)-G(0) (6.16) 
Jo 

+ P)fis)\\l + Ms)\\l + 1(1 - P,)f{s)\l^ + \r(s)\l + \\bg}ds. 

Prom (|6.10p . (|6.14p and (j6.16p . we prove the lemma by choosing e sufficiently small. ■ 

Proof of Proposition 16.11 The proof is parallel to the steady case, proof of Proposition 
13.11 We start with the approximating sequence, with = fo, 

dtf'-' + v-y.f+^ + uf+^-Kf = g, r+i(0) = /o, (6.17) 

J 



Step 1. We first show that for fixed j, f ^ P as ^ — )• oo. Notice that by the compactness of 
ti^fi ^ ^\\\ fe.+i\\2 I II ^£ii2i I ri \\\ fe.+i\\2 I w ft]^ ^ 
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We first take an inner product (using Green's identity) with f^^^ in ()6.17p . With the same 
choice of Cj in (j3.6p and using the boundary condition 



2 + (i-^) n\f^Hs)\\i+ f\r\s)\i^ds 











< e 



s)\\tds + 











+ C, f\r\l + Ce f U^^f \\f\s)\\lds+ f\\g{s)\\lds+^^^~^^' 
Jo Jo 1<*<^+1 Jo 



0||2- 



< e / \\f{s)\\tds+ il--r + - / \f{s)\i^ds 
Jo I J J i Jo 

+C, f\T\l + C, [\ ma^ W f{s)\\lds+ f\\g{s)\\lds+\\fo\\l 
Jo Jo 1<*<^+1 Jo 



< max 



1 - e 
t 



1^9 1 



l\\ns)\\lds+ f\f{s)\l^ds 
Jo 



+C, 1^ \r\i + Ce 



J J 

\ma^ \\r{s)\\lds+ f\\g{s)\\lds + 

1<«<^+1 Jo 



oiii- 



We choose e > sufficiently small, so that max|Y^, (1 — i)^ ~^ } — 1 ~ ^- Now we 
iterate again to have 



< (1 - ef {{l-e) f \\f-\s)\\lds + f \f{s)\l^ds 
Jo Jo 



t ft ft 

2 , „ II #U„Ml2 j„ I / ll„/„M|2, 



+ [{l-e) + l]{C, \r\i + Cej^ ^maxj|r(.)||^d. + + ||/o||^} 



< {l-eY^'[{l-e)[\\f{s)\\lds + j\fX, 

r ^ Tlrli + C, r max ||r(.)||i(i.+ f\\g{s)\\lds + ml\. 
e I Jo Jo i<*<^+i Jo 



(6.18) 



We therefore have, from = fo 



,maxj|r(t)||i 



{ t\r\l+ r max \\r{s)gds+ f 
^ Jo Jo i<*<^+i Jo 



{s)\\lds+\\fo\ 



+tmi+t\fo\i+}. 



By Gronwall's lemma we have 



max \\nt)\\i<e,,t 

l<i<£+l 



\r\l + 



{s)\\lds + \\fo\\l + t\\fo\\l + t\fo\lA, 
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where t > is fixed. This in turns leads to 



l<i<£+l 



max^^||r(t)||i+ / \\ns)\\l+ / \fis)\lds 



'0 Jo 

ft ft 



^''Alo ^'^'^l 11^(^)11''^" + + 
Upon taking the difference, we have 

dt[f+' -f]+v- V.[/+i - f] + u[f+' - f] = K[f - (6.19) 

with [f+^ - fm = and -/! = (!- i)P^[/ - f~']. Applying to - f 

yields 

\\f+\t) - f\t)h + r \\f^\s) - ns)\\l + f \f+\s) - f{s)\l_,ds 
Jo Jo 

< (1 - eY 1^ \f - f\l+ + ^ I jJI* max^ Wfis) - f-Hs)\\lds'j ■ 

From (j6.18p . this implies, from the Gronwall's lemma, that, by taking maximum over 1 < i < 
£ + 1, 



\\f^Ht)-f{t)h <e,,t {l-eyflf-f 

Jo 



t 

1 f0\2 

2,+ ' 



' \\f^\s) - f{s)\\l + f \f^\s) - f{s)\l^ds S.,t (1 - ef f - f\l^. 
Jo Jo 

Therefore, is a Cauchy sequence and — t- so that is the solution of 

dtf + v-V,f + Lf = g, /(0) = /o, /_ = (1 - + r. (6.20) 

Step 2. We let j — ?■ oo. Upon using Green's identity and the boundary condition and p.Sp . 
for any rj > 0, 

\\fml+ f'\\{i-p)f{s)\\ids+ f\{i-p,)f{s)\i^ds 

Jo Jo 

< 7?/ \P^fis)\l+ds + cJ' \r\l+ f Ms)\\lds + \\fo\\l (6.21) 
Jo Jo Jo 

Note that, from Ukai's trace theorem (Lemma 12. ip and ()6.20p 

f\P^f{s)\l+ds < f\f\s)1^4ds (6.22) 
Jo Jo 

< t\\fis)gds+ f\\d,[ff + vVM^]% 
Jo Jo 

'\\fis)gds+ f\{Lf,f)\+ fus)\\lds 
Jo Jo 

t rt rt 



< [ \\fis)gds+[ \\{I-P)fis)\\l+ [ Ms)\\lds. 
Jo Jo Jo 
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But from (j6.22p . choosing rj small, we obtain from (j6.2ip and (j6.22p that 



< / \\P{s)\\lds+ / \r\l+ \\g{s)\\lds + W\'^ 



10 JO 
It follows from the Gronwall's lemma that 



2- 



\Wml <t f\r\l + I ||5(^)||^ci.+ ||/o||i, (6.23) 







and hence 



\\pml+ f\\{i-p)f{s)\\ids+ f\{i-p,)f{s)\i^ds 

Jo Jo 

\r\l+ fUs)\\lds + \\M\l (6.24) 
Jo 

Since \\PP{s)\\l < ||/^'(s)||i, integrating (lU:^ from to t and combining with we have 

r \\Pf{s)\\lds + f \P^f{s)\2,+ds (6.25) 
Jo Jo 

<t f\\{I-P)f{s)tds+ [\\r{s)\l + \\g{s)g}ds+\\fo\\l 



Add (jOij) and ([OS]) together to have 



t 

fimi|2 + / II fi 



{s)\\lds + r \fis)\lds <t f \t\1 + f Us)\\lds + ml (6.26) 
Jo Jo Jo 



By taking a weak limit, we obtain a weak solution / to (|6.2p with the same bound (j6.26p . 
Taking difference, we have 

dt[P -f]+v - /] + L[P - /] = 0, [P - /]_ = P^[P - /] + -^P^f\ (6.27) 

with — /](0) = 0. Applying the same estimate (|6.26p with r = jP'yf'' we obtain 

\\P{t)-fml+ f\\P{s)-f{s)\\lds+ f\p{s)-f{s)\lds 
Jo Jo 

<t- f\P,P? <t-^o. 

J Jo J 
We thus construct / as a solution to (16. 2p . 

Remark 6.3 In both Step 1 and Step 2, we do not need the zero mass constraint which is 
instead essential for next step. 
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Step 3. Decay estimate. 

To conclude our proposition, let 

We multiply (|6.2p by e'^*, so that y satisfies 

dtV + v- V^y + Ly = \y + e^'g, y\^_ = P^y+ + e^V. (6.28) 
We apply Green's identity and (|6.22p . (|6.2ip together with ^ A to obtain 

\\y{t)\\l+ f\\{l-V)y{s)\\l+ f \{l - R^)y{s)\''+ (6.29) 
Jo Jo 

< A fus)g + \\ym\l + C, fe'y\l+ f e^^\\9{s)\\lds. 
Jo Jo Jo 

Prom (|6.2p we know that 

[[ yVJ^= [[ (Ay + e^'g)^ = 0, [ e^'r^dj = 0. 
JJnxR^ JJnxR^ J7_ 

Applying Lemma 16.21 to (|6.28p , we deduce 

t\\Py{s)\\lds<G{t)-G{0) (6.30) 

^0 

+ r\\{I-P)y{s)\\lds+ fe'^Wggds (6.31) 
Jo Jo 

+A / \\y\\lds+ [ {\{l-R^)y{s)\l+ + e^y\l}ds, 
Jo Jo 

where G{t) < ||y(i)ll2- Multiplying with a small constant x (|6.30p + (|6.29p to obtain, for some 

iUml - eG{t)] + e 11(1 - P)y{s)\\l + ||Py(.)||^} + " ^7)^(^)111+ 

< CX f \\y{s)\\l + [\\ym\l - eGm + G,,, f e^^\r\l+ f e^^\\g{s)\\lds. 
Jo Jo Jo 

By further choosing A <C e, and ||(I — P)2/(s)||^ + ||Py(s)||^ > ||y(s)||2 for hard potentials, we 
conclude our proposition. ■ 



7 L°° Stability and Non-Negativity 

To conclude the proof of Theorem 1 1.3 1 we need L°° estimates. They are provided, in the linear 
case, by next 
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Proposition 7.1 Let ||tt;p/o||oo + \{v)wpr\oo + ||tt;p5r||oo < +00 and JJ = J^r^ 
11 foy/J^ = 0- Then the solution f to i6.2\) satisfies 



fit)\\oo + \wpf{t)\oo <e ^'{ll'Wp/olloo +supe^^||tt;p5r||oo + / e^'\{v)wpr{s)\oods} 



Furthermore, «/ /o|7„ = -P7/0 + '"O; fo,^ o,nd g are continuous, then f{t,x,v) is continuous 
away from D. In particular, it Q is convex then 2) = 70. 

Proof. We use exactly the same approximating sequence and repeat step 1 to step 3 of 
Proposition 16. 1 1 to show that f^, f are bounded, via Lemma [4. 4[ We denote = Wgf^, where 
Wg is scaled weight in (14. Ij) . Rewrite (I6.17P as 



Kwgh +Wgg, 



(7.1) 



1 



"WgKV) Jn{x)-v'>0 



h {t, X, V )wg{v')da + w^r. 



Step l.We take £ — )• 00 in L°° . Upon integrating over the characteristic lines ^ = v, and 
^ = and the boundary condition repeatedly, along the stochastic cycles, we obtain (by 
replacing 1 — i with 1 and e = 0) that (|4.8|) is valid for Therefore, for £ > 2g, by 

Lemma 14.41 that 

sup e~^'^\h^~^^(s,x,v)\ 

0<s<To 

< i max sup {e'^^h'-^sn^j + Cik) max r \\f-\s)hds 

o 0<l<2k o<s<To l<l<2k Jq 

is)\oo} 

= - max sup |e 2 *||/i'^~'(s)||ool + C(A;) max / 

8 0<l<2ko<s<To i<l<2kjQ 



+ ||/io||oo + k 



sup <e 2 \Wgr[ 

0<s<To 



+ sup ■{ e"2 

0<s<To 



Wggis) 


00 } 


{v) 





'-Hs)hds 



+D, 



where Tq = g'^^^, g = k are chosen to be sufficiently large but fixed. Now this is valid for all 
£ > 2k. By induction on £, we can iterate such bound for £ + 2, ....£ + 2k to obtain 



sup e^'\\h^+\s) 

0<s<To 



(7.2) 



< ^ max snp {e'^'\\h'+'-'{s)\\oo} + C{k) max [ ' \\f-'{s)\\2ds + D 

< ^ max sup {e^"||/i^+^-'-^(s)||oo} + 2C(A;) max 1^° \\f-^(s)\\2ds + 2D 

4 l</<2fco<s<ro -2fe<i<2fc Jo 



< ^ max sup {e^-^'\\h^-^{s)\\^} + {i + l)C{k){ max / \\f-'{s)\\2ds + D}. 



4 l</<2feo<s<To 



-2k<l<2k 



To 



rl-l/ 
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Taking maximum over i = 1, ...2k, and by induction on 

max sup e^''||/i^+^~'(s)||oo 

l</<2fco<s<To 



< ^ max sup {e'^^\\h'-\s)\\oo} + C{k) 

4 l<i<2fco<s<To 



max 

-2fc<K2fc 



To 



\\f-\s)hds 



+ D 



1 



< max sup {e^^||/i^-i-'(s)||oo} + (l+ , 

42 l<«<2fco<s<To 4 



\){c{k)\ 



max 

'2k<l<2k 



To 



\\f-'{s)hds] +d] 



< 



max sup {e'i^\\h'{s)\\^} + {cik)\ max ^ \\f-\s)hds\ + d} , 

/^l/ZK l<Z<2fc o<g<jjj L l-2k<l<2k Jg J J 



where we may assume that ^ is a multiple of k. Now for maxi<;<2fc ||^'||ooi we can use (|4.8 
for A; = 1 repeatedly for /i^'^ — ). /i^^-i ...—)• /ig to obtain 



max sup {e"""!!/! (s)||oo} <fc ||/io||oo+ sup <e2''\wgr{s)\oo>+ sup <e 

l<'<2A:o<s<To 0<s<To 0<s<ro I 



2 * 



We therefore conclude that 



max sup e^"*!]/!^"^^ '(s) 
l<«<2fco<s<To 



feC(A;){ max [ ' \\f-^{s)hds + d} 

L — 2fc<(<2fc Jg J 

+ ||/io||oo+ sup |e^''|'u;gr(s)|oo[ + sup <^ 
o<s<To J o<s<ro L 



6 2* 



Now maxi<;<oo ||/'||2 is bounded by step 1 in the proof of Proposition 16. H || • II2 is bounded 
by \\wg ■ lloo for /3 > 3 and \r\2 is bounded by \wg{v)r{s)\oo. Hence, there is a limit (unique) 
solution hi ^ h = Wgf G Furthermore, h satisfies ()4.8p with h^'^^ = h. By subtracting 
h^~^^ — /i in (14.81) with D = 0, we obtain from Lemma 14.41 for h^^^ — h and 



max sup e 2 ''ll/i^'*"^ \s) - h(s) 
l<«<2fco<s<To 



< 



1772fc,™, sup {e 2''^||/i'(s) - /i(s)||oo} 



-{C(A.) 



max 

'2k<l<2k 



To 



\\f-^is)-f{s)hds]. 



From step 1 of Lemma WA] we deduce supg<g<2|j \\h^~^^ '(s) — /i(s)||oo — >■ for £ large. 



Step 2. We take j — )• cxd. Let to be the solution to (j6.20p and integrate along 
repeatedly. Again (|4.8p is valid for = (replacing (1 — i) by 1 and e = 0) and / = 0. 



dx 
dt 



dv 
dt 
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Lemma 14.41 implies 



sup e'^^\h^{s)\U 

0<s<To 
1 



<- sup {e^^\\h^{s)\\oo} + Cik) / \\fis)\\2ds 

° 0<s<To 



To 



+ 



+ k 



sup \e2''\wgr{s)\oo\ 
0<s<To ^ J 



+ sup <j e 2 ' 

o<s<ro 



(v) 



with To = g'^^^, Q = k, and g sufficiently large but fixed. Therefore, by an induction over j, 
sup e^^||/i^'(s)|U < Cik) r \\P{s)\\2ds + \\h4^ 

0<s<To JO 



+k 



sup \e'^''\wgr{s)\oo\ 
o<s<ro 



0<s<To 





^i'£)5(■5) 






{v) 


oo 



Since Jq^° ||P(s)||2(is is bounded from step 2 of Proposition [6TT] (with no mass constraint), this 
implies that ||/i-'||oo is uniformly bounded and we obtain a (unique) solution h = Wgf S L°°. 
Taking the difference, we have 



dt[y-h] + vV^[h^ -h] + u[h^ -h] = K^^[h^ -hi 



[h^ - h]_ 



Jn{xyv'>0 
1 1 



[h^ — h]{t,x,v')wg{v')da 



j Wq{v) Jn{xyv'>0 



y {t, X, v')wg{v')da. 



We regard j:;^^ /n(a;)-i)'>o ^' ^ Lemma So Lemma 

that 



implies 



sup e'^'Wh^is) - h{s)\\^ 

0<s<To 

<C(A;) r''||/^-(s)-/(s)||2<is + fc sup \w,r{s)\ 

Jo 0<s<To 



<i 

oo ^ . ) 



which goes to zero as j to oo. 

We now obtain a L°° solution /i = Wp/ to (|6.2p . Since we have L°° convergence at each 
step, we deduce that h is continuous away from D. 

To obtain decay estimate, we integrate (|6.2p over the trajectory. The estimate ()4.8p is 
valid for = h so that from Lemma 14.41 for g > sufficiently large we have: 



1^0 1^0 

sup e-''||/i(s)||oo + sup e-''\h{s)\c 

0<s<To 0<s<To 



(7.3) 



< 



+ C{k) 



sup 

0<s<To 



+ sup e 2 

0<s<To 
To 



(v) 



\\fm2ds 
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where Tq = g^^^, g = k and g sufficiently large. Letting s = Tq in (|7.3p . we obtain 



||/i(ro)||oo + |M?o)|oo <e-^^o||/io 



(7.4) 



+C{k) 



sup {\wgr{s)\oc} + sup 

0<s<To 0<s<To 



Wgg{s) 



{v) 



+ 1^" \\fis)\\2ds 



Let 



R= sup \e2''\{v)wgr{s)\oo\ + sup \e2''\\wgg{s)\\oc\ + \\wgfo\\oo- 



0<s<oo ^ ^ 0<s<oo 

Since /3 > 3/2, Proposition 16.11 implies that 

Il/Wll2<. e-^'R. 

Let m be an integer. We look at times t = mTo. By induction on m, we have from the 
definition of R and (I7.4p for A ^ z^O) 



||/i([m + l]ro)||oo + |/i([m + l]ro)|, 
< e-'?^«||/i(mro)||oo 



+C{g)e sup {\w gr{s + mTo)\ oo} + sup < 

.0<s<To 0<s<To L 



t(;g5(s + mTo) 



{v) 



(■To 



+ \\f{s + mn)\\2ds 
Jo 



< e-^^«||/i(mro)||oo + C(f5)e-^^"-^ x 



sup e~^2~ {\wgr{s + mTQ)\oo} + sup e" 



0<s<To 



0<s<t 
/■To 



'u;g5(s + mTo) 



{v) 



+ 



uXTq 



e—\\fis + mTo)\\2ds 





< e-^^« ||/i(mro)|U + C(£.)(l + To)e-^^«-'^i? 

< e-2^^»||/i([m - l]ro)|U + C{g){l + T^)e-^^^^~^R 
+C{g){l + ro)e-^^Oe-^^Oe-^^^i? 

= e-2?'^o||/i([m - l]ro)|U + C{g){l + ro)e-^^«-^i?[l + e 

< e-3f^«||/i([m-2]ro)||oo 

+C{g){l + ro)e-^^« 2^R[1 + e"^^" + e'^^To] 



Toi 



< e-^^o||/.o||oo + C(,)(l + ro)e-^^-^i?ye-^-^^o 



< C'fc^j.o^A.ge 2 R. 



Combining with (|7.3p for < ti < Tq, we deduce that for any t = mTo + ti, \\h(t)\\oo + 
l^(*)|oo ^k,uo,\,g e~*i?. We deduce our proposition with w = Wg ioi g > sufficiently large 
but fixed. A simple scaling of g concludes the proof of the proposition. ■ 
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By using the above L°° estimate we can conclude the proof of Theorem II. 3[ 
Proof of Theorem 11.31 We consider the following iteration sequence 



f. 



^+1 us- fJ- 



f^{n ■ v)dv. 



7+ 



Clearly we have 



nxR3 



{L^fj' + r{f,f)}^dxdv = 0, 



7+ 



fUi = 0. 



Recall w^{v) = (1 + £>2|t;|2)f gCM' in g^J. Note that for < C < i, 



< 6 sup \\e^^'wgf{s)\\oc + \ sup \\e^2Wgf{s) 

0<s<t lO<s<t 



Using O for < C < CT' 



e 2 Wg{v) 



7+ 



< 6 sup |e 2 r (s) 

0<s<t 



By Proposition 17.1^ we deduce 

sup \\e^Wgf+^{s)\\oo + sup \e^Wgf^'^{s) 



0<s<t 



0<s<t 



< \\Wgfo\\oo+S sup ||e 2 l(;g/^(s)||oo + (5 sup 16 2-^^/^(5) 
0<s<t 0<s<t 



+ <{ sup ||e 2 tt)^/ (s) 

,0<s<t 



For 6 small, there exists a eo (uniform in 5) such that, if the initial data satisfy (ll.20p . then 

sup \\e~Wgf^"^'^{s)\\oo + sup \e~Wf,f^^^{s)\oo <\\wgfo\\oo- 

0<s<t 0<s<t 



By taking difference / — / , we deduce that 

dtif^' -f]+v- V.[/+i - f] + L[f^' - f] 

= L^,f. [f - f-'] + nf - f-\f) + nf~\f - f-' 

with f^^^ — = initially. Repeating the same argument, we obtain 
sup \\e'^w,[f+^-f]{s)\U+ sup \e'^w,[f+^-f]{s)\^ 



0<s<t 



0<s<t 

As 



< [<5+ sup We-^ Wef{s)\\oo+ sup We^Wgf ^(s)||oo] sup \\e2Wg[f-f \s) 

0<s<t 0<s<t 0<s<t 
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This implies that f^~^^ is a Cauchy sequence. The uniqueness is standard. 
Positivity. To conclude the positivity of F^, we need to show 

Fs + ^f{t) > 0, 

if initially Fs + ^/Jlfo > 0. To this end, we first assume the cross section B is bounded and 
hence u is bounded. We note that previous approximating sequence is not well suited to 
show the positivity. Therefore, we need to design a different iterative sequence. We use the 
following one: 

dtF'+' + vV.F'+' + u{F')F'+' = Q,,,r.{F\F% 

Ff+i = US [ F^in{x) ■ v)dv, 

where F^ = Fs + y^/I/o and 

iy{F) = / dw* / dujB{v — v^,,uj)F{v^:), (7.5) 

with initial condition F^+^(0) = + ^/Jlfo > 0. Clearly, such an iteration preserves the non- 
negativity. We need to show F^ is convergent to conclude the non negativity of the (unique!) 
hmit F{t) > 0. Writing F^+i = Fs + we have 

dtf^' + V • + uiv)f+' - Kf (7.6) 

f-^' = P,f + ^ fVJ^inix) . v)dv. 

Taking inner product (Green's identity) with Z^"*"^, from Ukai's trace Theorem, and {Kf^, f^^^) ^ 
ll/^lli + ll/^^^'^lli (bounded cross section B) and the boundary condition we have 

\\f^\t)\\l+ f\\f^'\\l+ f\f^% 
Jo Jo 

ft rt 



Jo 

< \\fof2 + [l + C6'] f\P,f\l+ fwfWl (7.7) 
Jo Jo 



+ [5+\\wfU + \\wf\\U j \\f 

Jo 

By (|6.22p and the equation (|7.6p . assuming maxi<;<£ ||w/'||oo < +oo we have 

f\p.f\i < rii/-^i+[5+iK^-^iioo+iK^-^iiL] r 

Jo Jo Jo Jo 

< max / WfWlds. 



wfwl 

(7.8) 
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Splitting 1 + CJ^ = 1- 5 + {6 + C6'^) in ([721) and using (US]) we have 

\\f'-'\t)\\l+ f\\f^'\\l+ f\f^M^ 
Jo Jo 

< ll/olli + f \P,f\lds + max f \\f\\lds. 



Then we iterate to obtain 



f^+l|2 



11/^-^^111+ f\\f^'\\l+ f\f 

Jo Jo 

< [1 - sf' + [1 + (1 - «)l + , Il/'(«)ll2} 

Taking maximum over !</<£ + !, 

rnax ||/(t)i<||/0||2+ f \P,f\l+ \\f{s)gds. 
By Gronwall's lemma, from = /q, 



(7.9) 



max ||/'(t)||2 <t ll/olb + tlP^/ob, (7.10) 

l<l<n+l 



uniformly bounded. 

We now apply Lemma 14.41 with 

g = L^,,J + nf, f\ r = ^ f^dj, 

for some Tq large, with = fo, using the same arguments as in proof of (j7.4p . to get 

sup e 2 \wgj (s, X, v)\ 

0<s<Tq 

< ^ max sup {e^''||tt;p/+^~'(s)||oo} + ||^«e/o||oo 

4 l<«<2A:o<s<To 



-Tn 



(5 max sup je'a '(s)|oo| + max sup je^a" '^luj^Z+i '(s)|oo| 

l<'<^0<s<To J 0<«<2A;o<s<To 

To 

e-1. 



+C(ro) max / \\r-\s)hds. 



where i > 2k. For 

5+ max sup ie^"|?i;„/+^-'(s)|oo j < 1, (7.11) 
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small, we obtain from /? > 3, 

sup e^^\wQf^^{s,x,v)\ < I- max sup {e^^\\wgf^'^~^{s)\\^} + CTo\\wJo\\oo- 

0<s<t ^ \<l<2k o<s<i 

Hence we obtain 

max sup e'^^\wj'-^'^~\s,x,v)\<\ may. sup {e^*||tt;g/^+^"'(.s)||oo} + Cfclhe/o||oo, 
i<«<2fco<s<ro ^ i<'<2fco<s<ro 

and 

max sup e^''|u;g/+^"'(s,x,f)| < max sup {e^''||'W£,/(s)||oo} + 2Cfc||i(;g/o||oo- 

l<'<2fco<s<To 2v l<'<2fco<s<ro 

But for 1 < / < 2/c, we use Lemma with A; = 1 to get 

max sup {e 2 "^Wq^ (s)||oo} < C'fc||'«^£,/o||oo, (7.12) 

l<l<2k 0<s<To 

SO that (jT.lip is valid as long as ||tt;g/o||oo is sufficiently small. We therefore obtain uniform 
bound 

max sup \\wgf{s)\\oo < CfcHiy^./oHoo- 

l<'<^0<s<To 

This leads to Wgf^ Wgf G L°°. Furthermore, / satisfies (j7.6p with /^+^ = = f- 
Therefore f^~^^ — f satisfies 

{dt + vV, + i^ivmf-^' - /] - K[f -f] = R', - /](0) = 0, 

- /]+ = P,[f -f] + [ [f - f]^{n{x) . v)dv, 

where 

\iR',f^' - f)\ < {S + Wwj'Woo + \\w,f\\l + \\wj\\^ + \\wj\\l}{\\f+^ - (7.13) 
Notice that from (f71^ 

max II ij;^/ 1 1 00, ||w^£)/||oo < C'^||«^£)/o||oo. 
With the same proof of (|7.8p 

/V7[/'-/]l2< max l|r-/||id.. 
Combining the two above estimates we have 

\\f^Ht)-m\\i+ f\\f^'-f\\i+ f\f^'-f\i^ 

Jo Jo 
<[l-6] f \P,[f - f]\lds + max f \\f - f\\lds. 
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We iterate as for (j7.9p to obtain 

max writ) - f{t)\\l < r nmx Wf - fgds. 
i<i<(:+i Jq i<i<e+i 

Then applying Gronwall's lemma and choosing small > we conclude that is Cauchy 
in [0, t] and P — )• / by uniqueness and F{t) > for < t < to- We can repeat this process to 
show F{t) > 0, for ah t > 0. 

Finally, we take Fq ~ /i, so that Fq > and hence F(t) > 0. Therefore lim^^+oo I^{t) = 
Fg > 0. This completes the proof in the case B bounded. To remove this limitation, we use 
a cut-off procedure as before and reduce to the bounded case where previous result holds. 
Then we pass to the limit in the cut-off using the a priori bounds and uniqueness. ■ 
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